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A number of radiation shielding problems involve calculations of the response of an 
isotropic detector to radiation of arbitrary angular distribution from uniform finite plane 
sources. Series expansion methods previously used for treating the rectangular source are 
here extended to include the circular disk source with detector off-axis. These methods 
involve the family of integrals f 8 (cos 6 dS/r 2 )Pi (cos 0) and the integral J*s (dS/r 2 ) exp ( — fir) 
where is obliquity with respect to the normal to the disk surface S f r is the distance from 
an element of source area, dS, to the detector, and /x is the attenuation coefficient. Tabula- 
tions of the first type of integral facilitate use of Legendre expansion representations of 
radiation directional distributions. The second integral relates primarily to exponentially 
attenuated radiation from a plane isotropic disk source, but the expansion coefficient solution 
can be readily adapted to take into account a point isotropic buildup factor of the polynomial 

2V 
form ^j &,(jJir) 1 '. This adaptation applies equally well to the corresponding expansion 

t=o 
coefficients previously given for the rectangular source. Formulas and numerical results 
are presented. 



1. Introduction 

In radiation shielding, photometry and related 
studies, a frequently occurring problem is that of 
calculating the radiation field resulting from a finite 
plane source of radiation with arbitrary polar angular 
distribution g(6). A number of authors have in- 
vestigated the response of surface-type [l] 2 detectors 
to rectangular and circular Lambert's-law radiators. 
Smith and Storm [2] have derived a generalized 
series expansion method for calculating the response 
of an azimuthally symmetric detector off-axis from 
a disk source witli arbitrary surface angular dis- 
tribution, in which evaluation of the series terms 
involves successive differentiation of the detector 
response to a differential element of the source 
area. The resulting formulas usually consist of 
triple or quadruple summations in which one sum 
can be obtained by recursion methods. 

In a previous paper [3], a Legendre expansion 
method suggested by Berger and Lamkin [4] was 
applied to a rectangular source as part of a general 
study of the penetration of radiation into struc- 
tures [5]. In the Berger-Lamkin method, the geo- 
metric relationship of detector to source is character- 
ized by a set of Legendre coefficients, p h over which 
corresponding penetration-dependent coefficients, g h 
of a Legendre expansion of g(6) may be summed. 
This has some advantage over the results in reference 
[2] in that not only are the geometrical and penetra- 
tion parameters separated for economy of data 
tabulation, but use of the present Spencer-Fano [6] 
moments method provides most directly the g t 
coefficients, rather than the actual angular dis- 
tributions g{6). 



For small values of fit, where ju is the attenuation 
coefficient and t is the barrier thickness, the Legendre 
series solution does not converge adequately for the 
particular angular distribution for the unscattered 
radiation from a plane isotropic source 3 

0(0) = <AO) = W4tt) exp (— fit/cos 0)/cos 6, (1) 

in which a/iw is the source strength in units of 
radiation (e.g., photons) emitted per unit time from 
unit surface area into unit solid angle. The reason 
for this poor convergence is that, as jit-^O, g°(0)-+ 
(cr/47r)/cos 0, and a function of the form 1/cos 6 [is 
not represent able by the positive powers of cos 6 
appearing in the Legendre polynomials. An alterna- 
tive series expansion in powers of nt with geometry 
dependent coefficients, q n , was provided for calcu- 
lating the unscattered radiation component when 

0<m*<1. 

The present work extends the two expansion 
methods in reference [3] to include the case of a 
circular disk source with detector off-axis. A typical 
example of this situation would be the estimation 
of the neutron flux around the face of a cylindrical 
thermal column of a research reactor [7]. In addition, 
it will be shown that if the scattering properties of 
the medium are given as a point isotropic buildup 
factor [6] in polynomial form [8, 9] 



5(/*r)^SA(/^)S 
*=o 



(2) 
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rather than as angular harmonics, g h at the exit 
surface of the barrier, the q n expansion coefficients 
here and in reference [3] can often be used to evaluate 
the scattered as well as unscattered contribution to 
the radiation flux. 



3 Here and in what follows, g° and 1° will be used to denote the angular distribu- 
tion and intensity of direct radiation, i.e., radiation which arrives at the detector 
from the source after having undergone no scattering interactions with the in- 
tervening medium. 
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2. Arbitrary Angular Distribution; Legendre 
Expansion of the Angular Response Function 

2.1. Isotropic Detector (Flux) 

If the emission from each source element, dS, 
of the disk, is independent of azimuth and uniform 
over the source surface S, the radiation intensity at 
any point in space, as measured by a small spherically 
symmetric detector, may be written as the integral 



I=Sg(o).Mo)d(cQae). 



(3) 



Here 6 is the polar angle as indicated in figure 1. 
The angular distribution function gift) of radiation 
crossing the source -plane will now be defined as the 
radiation flux 4 (e.g., photons) traveling in direction 
6 per unit solid angle, unit time, and unit area normal 
to direction 0. $s(0), called the "angular response 
function" in reference [4], is the azimuthal angle 
subtended by the source area from the detector-to- 
source -plane vertical for a given 0-cone around this 
vertical. 

There are a variety of ways in which the integral 
(3) may be evaluated, depending on the nature of 
g{6) and S. The Legendre expansion method, es- 
pecially advantageous if g z coefficients happen to be 
already available, has been discussed in detail in 
references [3] and [4], but for convenience a brief 
summary is given as follows: 

The two functions, g(6) and $s(0), in (3) which are 
parametrically independent except in 0, may (if free 
of singularities) be expanded in Legendre harmonics 



as 



and 



1=1) 



-l)flf,-P,(costf) 



1=0 



(4) 



(5) 



* In this notation the radiation flux in direction d per unit solid angle, unit 
time and unit area of the source would be #(0)cos 0, i.e., g(d) may be identified with 
the angular distribution functions g(z,d) and <Mm) of references [4] and [7], re- 
spectively, but differs from (?(/?) of reference [2] according to #(0) = (?(/?) /cos 0. 



PLANE OF S 




DISK SOURCE, S 



Figure 1. Disk source and off-axis detector, showing polar 
angle and "angular response function" \p 8 (0). 



When this is done, the integral (3) may be written 
as the series 



I=T,(l+\)9i-Pi 



1=0 



(6) 



which converges rapidly for many cases of interest. 
The pi coefficients may be evaluated in the usual 
way according to 



Pl= J P,(cos 0) • ifc(0)d(eos 0), 



(7) 



or since ife(0)<i(cos 0) is a differential solid angle 

(8) 



Pi= J Pi (cos d)dtt s (d) 



As in reference [3], the \ps(0) and related p/s are 
dependent on the relative, rather than on the abso- 
lute dimensions, so that we may select one dimension, 
here the disk radius, as the unit of length. The 
quantities p and h, as indicated in figure 2, are then 
the distances off axis and off disk-plane, respectively, 
measured in disk radii. 

Since the polynomials P z (cos 6) are linear combina- 
tions of even or odd powers of (cos 6) it is convenient 
to break up (8) into components Fj(p,h), each over a 
particular power of cos 0, i.e. 



F j (p,h)= [ cos* edQ(6,p,h), 



(9) 



from which we reconstruct the desired functions 
Pi(pjh) according to 



Pi(p,h) = ^2a>ijFj(p,h), 



(10) 



where a t j are the coefficients of x j in the Legendre 
polynomial 

of order I. 

The differential solid angle, dQ(d,p,h), can be 
written in cylindrical coordinates h, R, and </> around 




DISK SOURCE ,S 



Figuke 2. Source- detector geometry, showing relevant param- 
eters p, h, 6, </>, r and R, in which linear dimensions are 
measured in radii of the source disk. 
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the perpendicular from detector to source-plane as 
dfl(0, P ,A)=cos 9 dS/r 2 =(h 2 +R 2 )' m hR dR </</>, (11) 

since cos 6=h/^h 2 +R 2 , dS=R dRcty, and r 2 =h 2 +R 2 . 
The integral in (9) may now be written explicitly as 
the double integral 5 



Fj(p,h) 

fl = p-lj0=-l 



,p,h) 



/2t{6,p,h) 



(h 2 +R 2 )- u+ * )/2 RdRdct>, 



(12) 



in which the limits ±^(0,p,A) of the azimuthal 
variable are indicated in figures 1 and 2. The 
angular response function \f/(6,p,h) is the quantity 
fs(0) defined following eq (3), in which the general- 
ized source-detector geometry symbol S has been 
replaced by the parameters p and h which character- 
ize the disk source situation. 



(13) 



Now, since 

^(d, P M)=2 co^ l {(p 2 +R 2 -l)/2 P R} , 
we may write (12) as the single integral 

Fj( P ,h)=h j+1 ( 2RdR(K 2 +R 2 )-«+ z)l2 

Xcos- 1 {(p 2 +7? 2 -l)/2 P K}. (14) 
Integration by parts and use of the substitution 
5 2 =p 2 +l-2p cos t (15) 

transforms (14) to the form 



{ 



'"-"I 



dt 



H i 
(p 2 -^ -2p cos t)(h 2 +p 2 + l-2 P cos 5 

dt *\ 



•r 



i k (16) 
(A 2 + p 2 + l-2pcos^) 2 J 



The first of the two integrals in (16) may be further 
reduced by use of the identity 



1 



a 2 (h 2 +a 2 ) j a 2 h 2j 






1 



1 



h 2j-2i (f l 2 +a 2 ) i 



(17) 



in which a 2 =p 2 +l — 2p cos t. 

Since at this point the evaluation of the integral 
for the case of oddj (and thus odd I) is quite different 
from that for even j, these cases will be (rented 
separately in subsections A and B, as follows: 



6 The limits of integration p— 1, p-(-l in (12) refer to the case l<p<°°. For 
0<p<l the integral should be split into the ranges 0<JR<1— p and 1— p</?<l+p. 
However, the resulting expressions arc identical. 



A. Even j, I 

Combining (16) and (17) we have for the F- 
integrals for even j } denoted as F v , where j=0, 

1,2..., 

F 2j ( P ,h)=r 



2j+l 



X 



r^-ir 

I A" J 



dt 



V 



( P 2 +l-2pcosOV^ 2 +P 2 +l-2p cos t 

i)i"§— r 



■j; 



s t) i+z/2 j 
(h 2 +p 2 + -2 P eos*)' +I7 ~ 2 J ^ 



o (A 2 + P 2 +l-2p cos 
dt 



The first, integral in (18) can be evaluated in terms 
of complete and incomplete elliptic integrals of the 
first and second kinds. However, in combining 
these F 2 /s according to (10) to give p 2 i, it can be 
shown that for Z>0 this integral cancels out, i.e., 
the polynomial 2a^(2j+l) _1 has the properties 



I>,(2j + l)-i= 



1 forZ=0 

_0 for Z>0, 1 even_ 



(19) 



The zeroth coefficient p (p,h) (which is the same 
as F (p,h)) is the response of an isotropic detector, 
at an off-axis position p,h, to radiation from a 
"perfectly diffuse" circular disk source, i.e., g(6) = 
constant. Fortunately, Po(p,h) is also the solid 
angle subtended by a circular disk from the point 
p,h, and in the latter role has been discussed and 
evaluated by a number of authors [10, 11]. If we 
use one of these tabulations 6 for p , we need only 
evaluate the integral appearing in (18) in the sum- 
mation and also in the final term. Factoring out 
the quantity (A 2 +p 2 +l) _(n+1/2) from this integral 
reduces it to the form 



s. 



-(*)=// 



(1-F cos t)- (n+m) dt, 



(20) 



in which n corresponds to the i and j in (18) and 
& 2 =2p/(A 2 +p 2 + 1) < 1. This family of integrals may 
be evaluated in terms of complete elliptic integrals 
of the first and second kinds [12]. 



S n (k) = 



S 1 (k)= 
and for n> 2, 



Vi+& 2 

9 



K\ 



(Viw 

( im\ 

(l-P)Vl+& 2 \\l+k 2 / 



(21) 
(22) 



(23) 



• For tabulations of po see, for example, references [10] and [11]. The param- 
eters p and h in this paper were explicitly selected in order to use directly the 
extensive table in reference [11]. 
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Extensive tabulations [13] of K and E exist in addi- 
tion to an excellent polynomial approximation by 
Hastings [14] suitable for use on a high-speed 
computer. 

Substituting (20) into (18) we have 

F 2 j(p, h) = {cancellation terms in p 2 i{p,h) for />0} 

h 2 



2j 



+i\W+p 2 +i) Sj(k) 






n 2_l ; / £2 \i + l/2 ^ 

The even-Z response function coefficients p 2 i(p,h) may 
then be built up from the F 2j (p,hYs in (24) as indi- 
cated in (10). Equation (10) may be written more 
explicitly for even I as 



j=0 



■iy/21-j 



rj^t-D'- 



j(p,h), 



(25) 



in which the quantity 

\r) = rl(n-r)l (26) 

coefficient of an nth degree binomial 

B. Oddj,l 

Denoting the odd-j ^-integrals in (16) as F 2j+l (p,h) 
where j=0, 1, 2, . . ., we have, using (17), 



is the rth 
expansion. 



F*.i+i(p,h) 



1 



2j+2 



{"-"IT 



dt 



(p 2 + 1- 
dt 



(h 2 +p 2 +l-2 P cos t) 



2p cos t) 



dt 



(h 2 +p 2 + l-2 P cos tV 



»}■ 



(27) 



Since the integrand denominators have integers, 
rather than half-integers, as exponents, closed-form 
algebraic expressions result, rather than quantities 
in terms of tabulated higher functions or infinite 
series expansions. Use of standard integral tables 
[15] to evaluate (27) gives 



F 2j +i(p 9 h)-- 



~2j+2 
h 2 



{ 



i 



W(A 2 -P 2 + l) 2 +4p%V ' \V(A 2 -p 2 +l) 2 +4pW 



-P 2 + 1) 2 +4 P 2 A 2 

— ±c 

A 2 U\. 



p-—l 



Y+i 

7 

/_ A 2 +p 2 +1 \\ 
, W(A 2 -p 2 +i) 2 +4pWi' 



-p 2 + l) 2 +4p 2 /i, 2 

i 2 +p 2 + l 



-p 2 + l) 2 +4p 2 A 2 , 
in which (cos 0) m is replaced by 

{(/> 2 +P 2 +D/V(A 2 -p 2 +i) 2 +4p% 2 r 



(28) 



in the usual Legem dre polynomials P n (cos 6) to give 
the quantities P f and P j in (28). 

Substituting (28) into (10), we have for the first 
two odd-l coefficients 



and 
Ps(p,h)-- 



Pl(p,h) 


<>- 


A 2 +p 2 -1 


V(/t 2 -p 2 +l) 2 +4p% 2 / 


=I{- 


1 1 


6A, 2 4-p 2 -l 


1+ vW 


-p 2 + l) 2 +4p 2 4 2 




5h i (h i + P 2 -l)(h i +p i + l) 



(29) 



[(A'-p* + l)*+4p*A*] 3 



*Y 



(30) 



The higher odd-Z terms can be written in a form com- 
parable to (25) as 



p 2 i+i(p,h)-- 






(-1)' 



/2Z-j + l\/4/-2j+2\ 
V j J\2l-j+l) 



F 2l - 2j+1 (p,h) (31) 



using the F 2j +i(p,h)'s from (28). 

The pi(p } h) coefficient in (29) is equivalent to the 
response of a plane detector parallel to the plane of a 
perfectly diffuse (Lambert's law) uniform circular 
disk source, and has been derived and discussed by a 
number of authors [16]. 

The F 2j (p,h) and F 2j+1 (p,h) expressions in (24) and 
(28), and the expressions in (25) and (31) for their 
respective Legendre combinations were programed 
for evaluation of the pi(p,h)'s on an IBM 704 
computer. The resulting numerical values for the 
Pi(p,h) are presented in table 1 for 0<p<10, 0.1 
<A<10 and 1 <£<13. The entries for Po(p,h) were 
taken from reference [11]. For p or h greater than 
10, the source may be treated for all practical pur- 
poses as a point anisotropic source with the same 
angular distribution, g(9) cos 0, as each surface element 
of the source. The Pi(p,h) coefficients in table 1 
may now be combined with g x coefficients as in 
equation (6) to give the radiation flux intercepted by 
a small spherical detector at coordinates p,h. 

As a check on these pi(p,h)'s, they may be inserted 
in (5) to reconstruct the angular response function 
\p(6,p,h). Three detector positions (p=0.5, A=0.5; 
P=l, h=l; p=2,h=l) are indicated in figure 3, and 
the exact angular response functions for these positions 
are given by the solid lines in figure 4. Summing 
over the first 10 terms in (5) gives us the approxi- 
mations shown by the dashed lines in figure 4. 

2.2. Cosine Detector 

As has been discussed in reference [3], a similar 
analysis of a finite uniform plane source geometry 
may be performed in the case of a surface-element 
detector with response proportional to the cosine 
of the angle of incidence of the radiation. For the 
case of a surface-element detector parallel to the 
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Table 1.* Coefficients pi(p,h) of the Legendre expansion (5) of the angular response function \f/ (9, p,h) 

Pi(p, h) 



/i = 0.1 


I \ 
\ 





0.2 


0.5 


0.8 


1.0 







5.6580 (0) 


5.6389 (0) 


5.5092 (0) 


4.9511 (0) 


2.7036 (0) 




1 


3.1105 (0) 


:<. 1079 (0) 


3.0872 (0) 


2.9346 (0) 


1.4924 (0) 




2 


3.0951 (-1) 


3. L855 (-1) 


3.7857 (-1) 


5.9042 (-1) 


1.6880 (-1) 




3 


-7.3913 (-1) 


-7.3534 (-1) 


-7.0554 (-1) 


-5.1185 (-1) 


-3.2419 (-1) 




4 


-2.2677 (-1) 


-2.3270 (-1) 


-2.6957 (-1) 


-3.3127 (-1) 


-9.7153 (-2) 




5 


3.3572 (-1) 


3.3122 (-1) 


2.9669 (-1) 


1.2095 (-1) 


1.4395 (-1) 




6 


1.8207 (-1) 


1.8596 (-1) 


2.0681 (-1) 


1.7067 (-1) 


6.9956 (-2) 




7 


-1.8038 (-1) 


-1.7550 (-1) 


-1.3959 (-1) 


-1.8423 (-2) 


-7.7337 (-2) 




8 


-1.5130 (-1) 


-1.5352 (-1) 


-1.6127 (-1) 


-7.5020 (-2) 


-5.2967 (-2) 




9 


1.0068 (-1) 


9.5697 (-2) 


6.1222 (-2) 


5.8812 (-3) 


4.6161 (-2) 




10 


1.2734 (-1) 


1.2811 (-1) 


1.2488 (-1) 


3.1191 (-2) 


4.2324 (-2) 




11 


-5.3333 (-2) 


-4.8521 (-2) 


-1.7867 (-2) 


-1.7343 (-2) 


-2.8446 (-2) 




12 


-1.0734 (-1) 


-1.0684 (-1) 


-9.4741 (-2) 


-2.2105 (-2) 


-3.4397 (-2) 




13 


2.2831 (-2) 


1.8417 (-2) 


-6.5421 (-3) 


2.4535 (-2) 


1.7852 (-2) 




\ P 

I \ 

\ 


1.2 


1.5 


2.0 


5.0 


10.0 







6.0618 (-1) 


1.7518 (-1) 


5.3816 (-2) 


2.6296 (-3) 


3.1768 (-4) 




1 


1.3539 (-1) 


1.9422 (-2) 


3.4561 (-3) 


5.4493 (-5) 


3.2047 (-6) 




o 


-2.4825 (-1) 


-8. 3956 (-2) 


-2.6553 (-2) 


-1.3131 (-3) 


-1.5879 (-4) 




3 


-1.7525 (-1) 


-2.8306 (-2) 


-5.1412 (-3) 


-8.1680 (-5) 


-4.8063 (-6) 




4 


1.0623 (-1) 


5 6815 (-2) 


1.9299 (-2) 


9.8181 (-4) 


1.1901 (-4) 




5 


L6643 (-1) 


3.3581 (-2) 


6.3309 (-3) 


1.0196 (-4) 


6.0060 (-6) 




6 


-1.4981 (-2) 


-3.9768 (-2) 


-1.5289 (-2) 


-8.1427 (-4) 


-9.9063 (-5) 




7 


-1.2740 (-1) 


-3.6289 (-2) 


-7.2270 (-3) 


-1.1873 (-4) 


-7.0039 (-6) 




8 


-3.5553 (-2) 


2.6539 (-2) 


1.2454 (-2) 


7.0782 (-4) 


8.6547 (-5) 




9 


7.8309 (-2) 


3.6858 (-2) 


7.9<X)6 (-3) 


1.3324 (-4) 


7.8748 (-6) 




10 


5.1249 (-2) 


-1.5553 (-2) 


-1.0188 (-2) 


-6.3174 (-4) 


-7.7741 (-5) 




11 


-3.6546 (-2) 


-3.5609 (-2) 


-8.3855 (-3) 


-1.4611 (-4) 


-8.6560 (-6) 




12 


-4.3885 (-2) 


6.3576 (-3) 


8.2513 (-3) 


5.7324 (-4) 


7.1095 (-5) 




13 


1.1524 (-2) 


3.2871 (-2) 


8.7016 (-3) 


1.5770 (-4) 


9.3693 (-6) 



h = 0.2 


\ 

\ P 
I \ 
\ 





0.2 


0.5 


0.8 


1.0 







5.0510 (0) 


5.0156 (0) 


4.7828 (0) 


3.9593 (0) 


2.4052 (0) 




1 


3.0208 (0) 


3.0113 (0) 


2.9394 (0) 


2.5521 (0) 


1.4145 (0) 




2 


5.9242 (-1) 


6.0653 (-1) 


6.9042 (-1) 


8.1628 (-1) 


2 6714 (-1) 




3 


-6.0996 (-1) 


-5.9725 (-1) 


-5.0554 (-1) 


-1.7057 (-1) 


-2.5690 (-1) 




4 


-4.0444 (-1) 


1.0939 (-1) 


-4.2399 (-1) 


-2.4971 (-1) 


-1.4043 (-1) 




5 


L.8601 (-1) 


1.7310 (-1) 


9.0334 (-2) 


-7.4080 (-3) 


9.3690 (-2) 




C) 


2.8843 (-1) 


2.8648 (-1) 


2.5294 (-1) 


5.7523 (-2) 


9.3989 (-2) 




7 


-2.4333 (-2) 


-1.3696 (-2) 


4.0103 (-2) 


-3.7853 (-2) 


-3.5202 (-2) 




8 


-2.0024 (-1) 


-1.9326 (-1) 


-1.3335 (-1) 


-5.9687 (-2) 


-6.4854 (-2) 




9 


-4.9727 (-2) 


-5.6315 (-2) 


-7.1927 (-2) 


2.6729 (-2) 


1.0905 (-2) 




10 


1.2879 (-1) 


1.1895 (-1) 


5.5438 (-2) 


7.2910 (-2) 


I 6952 (-2) 




11 


8.1495 (-2) 


N. 3107 (-2) 


6.2859 (-2) 


1.8096 (-2) 


1.7315 (-3) 




12 


-7.0789 (-2) 


-6.0437 (-2) 


-1.2489 (-2) 


-4.1686 (-2) 


-3.3718 (-2) 




13 


-8.8823 (-2) 


-8.5486 (-2) 


-4.0990 (-2) 


-2.8918 (-2) 


-7.7467 (-3) 




\ 
\ P 

l \ 




1.2 


1.5 


2.0 


5.0 


10.0 




9.7005 (-1) 


3.3005 (-1) 


1.0554 (-1) 


5.2489 (-3) 


6.3507 (-4) 




1 


3.6408 (-1) 


7.0472 (-2) 


1.3424 (-2) 


2.1740 (-4) 


1.2811 (-5) 




2 


-2.4386 (-1) 


-1.3978 (-1) 


-5.0043 (-2) 


-2.6108 (-3) 


-3.1715 (-4) 




3 


-3.4860 (-1) 


-9.4767 (-2) 


-1.9485 (-2) 


-3.2513 (-4) 


-1.9203 (-5) 




4 


-5.6872 (-2) 


6.5955 (-2) 


3.2928 (-2) 


1.9343 (-3) 


2.3718 (-4) 




5 


1.7291 (-1) 


9.6607 (-2) 


2.2940 (-2) 


4.0426 (-4) 


2.3975 (-5) 




6 


1.2253 (-1) 


-1.5020 (-2) 


-2.1798 (-2) 


-1.5809 (-3) 


-1.9676 (-4) 




7 


-3.1183 (-2) 


-8.2335 (-2) 


-2.4499 (-2) 


-4.6802 (-4) 


-2.7921 (-5) 




8 


-6.7218 (-2) 


-1.9212 (-2) 


1.2973 (-2) 


1.3466 (-3) 


1.7111 (-4) 




9 


-2.8126 (-3) 


5.8673 (-2) 


2.4463 (-2) 


5.2123 (-4) 


3.1338 (-5) 




10 


2.3823 (-2) 


3.7369 (-2) 


-5.5921 (-3) 


-1.1705 (-3) 


-1.5279 (-4) 




11 


-1.2139 (-2) 


-3.2733 (-2) 


-2.3072 (-2) 


-5.6621 (-4) 


-3.4373 (-5) 




12 


-2.6654 (-2) 


-4.1470 (-2) 


-5.4447 (-4) 


1.0276 (-3) 


1.3873 (-4) 




13 


9.2864 (-3) 


1.0522 (-2) 


2.0585 (-2) 


6.0426 (-4) 


3.7109 (-5) 



/i=0.5 


\ 

\ P 
I \ 

\ 





0.2 


0.5 


0.8 


1.0 







3. 4733 


(0) 


3.4184 (0) 


3.1022 (0) 


2.4123 (0) 


1.7687 (0) 




1 


2. 5133 


(0) 


2.4802 (0) 


2.2733 (0) 


1.7426 (0) 


1.1898 (0) 




2 


1.1240 


(0) 


1.1232 (0) 


1.0865 (0) 


8.2414 (-1) 


4.2719 (-1) 




3 


0. 0000 




2.3936 (-2) 


1.3416 (-1) 


1.5252 (-1) 


-7.3353 (-2) 




4 


-4. 4959 


(-1) 


-4.2364 (-1) 


-2.6953 (-1) 


-6.8444 (-2) 


-1.6341 (-1) 




5 


-3. 0159 


(-1) 


-2.9631 (-1) 


-2.1987 (-1) 


6.2302 (-3) 


-2.6599 (-2) 







4. 4959 


(-2) 


2.5241 (-2) 


-2.9195 (-2) 


9.8306 (-2) 


7.5291 (-2) 




7 


2.2117 


(-D 


1.9522 (-1) 


6.6073 (-2) 


7.1346 (-2) 


4.7875 (-2) 




8 


1.3488 


(-D 


1.2808 (-1) 


3.9757 (-2) 


-1.6869 (-2) 


-2.4193 (-2) 




9 


-5. 2276 


(-2) 


-3.2603 (-2) 


-8.3723 (-3) 


-5.7460 (-2) 


-3.9902 (-2) 




10 


-1.3847 


(-D 


-1. 1162 (-1) 


-8.4522 (-3) 


-2.8599 (-2) 


-1.1920 (-4) 




11 


-6.9165 


(-2) 


-6.2655 (-2) 


1.9010 (-2) 


8.6621 (-3) 


2.8202 (-2) 




12 


5. 179:< 


(-2) 


3.0551 (-2) 


2. 1923 (-2) 


5.9185 (-3) 


1.2729 (-2) 




13 


9. 5706 


(-2) 


6.8208 (-2) 


-7.9961 (-3) 


-2.0240 (-2) 


-1.4838 (-2) 
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Pi(p, h) 



h=0.5 


\ 
\ p 
I \ 

\ 


1.2 


1.5 


2.0 


5.0 


10.0 







1.1661 (0) 


5.9692 (-1) 


2.3242 (-1) 


1.2945 (-2) 


1.5826 (-3) 




1 


6.7870 (-1) 


2.6382 (-1) 


6.9462 (-2) 


1.3339 (-3) 


7.9725 (-5) 




2 


6.7922 (-2) 


-1.0857 (-1) 


-8.3311 (-2) 


-6.2644 (-3) 


-7.8526 (-4) 




3 


-2.7309 (-1) 


-2.3352 (-1) 


-8.6034 (-2) 


-1.9643 (-3) 


-1.1908 (-4) 




4 


-2.4935 (-1) 


-9.7788 (-2) 


1.5867 (-2) 


4.3407 (-3) 


5.7842 (-4) 




5 


-6.2514 (-2) 


8.1473 (-2) 


7.3710 (-2) 


2.3747 (-3) 


1.4771 (-4) 




6 


4.6779 (-2) 


1.1907 (-1) 


2.6715 (-2) 


-3.1647 (-3) 


-4.6834 (-4) 




7 


2.2690 (-2) 


3.0262 (-2) 


-4.4313 (-2) 


-2.6379 (-3) 


-1.7041 (-4) 




8 


-3.1217 (-2) 


-4.9324 (-2) 


-4.5078 (-2) 


2.2555 (-3) 


3.9367 (-4) 




9 


-2.5263 (-2) 


-4.4582 (-2) 


1.1501 (-2) 


2.7791 (-3) 


1.8892 (-4) 




10 


2.1162 (-2) 


9.5174 (-4) 


4.2228 (-2) 


-1.4823 (-3) 


-3.3618 (-4) 




11 


4.0843 (-2) 


1.7678 (-2) 


1.3162 (-2) 


-2.8121 (-3) 


-2.0405 (-4) 




12 


1.7192 (-2) 


6.0957 (-4) 


-2.6483 (-2) 


8.0073 (-4) 


2.8842 (-4) 




13 


-8.5251 (-3) 


-9.0943 (-3) 


-2.4048 (-2) 


2.7473 (-3) 


2.1626 (-4) 



ft = 1.0 


\ 

\ 





0.2 


0.5 


0.8 


1.0 







1.8403 (0) 


1.8071 (0) 


1.6371 (0) 


1.3488 (0) 


1.1227 (0) 




1 


1.5708 (0) 


1.5394 (0) 


1.3760 (0) 


1.0921 (0) 


8.6832 (-1) 




2 


1.1107 (0) 


1.0855 (0) 


9.4773 (-1) 


6.9040 (-1) 


4.8242 (-1) 




3 


5.8905 (-1) 


5.7669 (-1) 


4.9529 (-1) 


3.0375 (-1) 


1.3416 (-1) 




4 


1.3884 (-1) 


1.4450 (-1) 


1.4522 (-1) 


5.4083 (-2) 


-5.9580 (-2) 




5 


-1.4726 (-1) 


-1.2475 (-1) 


-4.2344 (-2) 


-2.8548 (-2) 


-8.8156 (-2) 




6 


-2.4297 (-1) 


-2.1313 (-1) 


-8.4731 (-2) 


5.7923 (-4) 


-2.6783 (-2) 




7 


-1.8715 (-1) 


-1.6474 (-1) 


-4.8632 (-2) 


5.4400 (-2) 


3.2267 (-2) 




8 


-5.8573 (-2) 


-5.6131 (-2) 


-4.4288 (-3) 


7.2684 (-2) 


4.1359 (-2) 




9 


6.0592 (-2) 


4.0296 (-2) 


9.2186 (-3) 


4.7889 (-2) 


1.1747 (-2) 




10 


1.1660 (-1) 


8.3307 (-2) 


-6.4612 (-3) 


9.6586 (-3) 


-1.7466 (-2) 




11 


9.8846 (-2) 


7.0358 (-2) 


-2.7229 (-2) 


-1.0687 (-2) 


-2.0649 (-2) 




12 


3.3964 (-2) 


2.6471 (-2) 


-3.1179 (-2) 


-5.3664 (-3) 


-2.9665 (-3) 




13 


-3.4778 (-2) 


-1.6224 (-2) 


-1.4585 (-2) 


1.0427 (-2) 


1.3218 (-2) 




v 

I \ 
\ 


1.2 


1.5 


2.0 


5.0 


10.0 







9.0037 (-1) 


6.1910 (-1) 


3.2580 (-1) 


2.4692 (-2) 


3.1293 (-3) 




i 


6.5297 (-1) 


3.9677 (-1) 


1.6583 (-1) 


5.0026 (-3) 


3.1407 (-4) 




2 


2.8932 (-1) 


8.8340 (-2) 


-3.1552 (-2) 


-1.0811 (-2) 


-1.5173 (-3) 




3 


-1.6324 (-2) 


-1.3518 (-1) 


-1.2927 (-1) 


-6.9761 (-3) 


-4.6313 (-4) 




4 


-1.5564 (-1) 


-1.8651 (-1) 


-8.8713 (-2) 


5.6155 (-3) 


1.0564 (-3) 




5 


-1.3742 (-1) 


-1.0432 (-1) 


1.4583 (-2) 


7.6054 (-3) 


5.6125 (-4) 




6 


-5.2771 (-2) 


3.1764 (-3) 


7.5614 (-2) 


-1.8398 (-3) 


-7.7684 (-4) 




7 


6.4034 (-3) 


5.2785 (-2) 


5.4384 (-2) 


-7.1757 (-3) 


-6.2568 (-4) 




8 


7.7352 (-3) 


3.7389 (-2) 


-6.5513 (-3) 


-1.1179 (-3) 


5.6176 (-4) 




9 


-2.1418 (-2) 


3.4071 (-3) 


-4.2034 (-2) 


5.9087 (-3) 


6.6240 (-4) 




10 


-3.8294 (-2) 


-7.6326 (-3) 


-2.9432 (-2) 


3.2646 (-3) 


-3.7906 (-4) 




11 


-2.6571 (-2) 


6.0364 (-3) 


3.1900 (-3) 


-4.0690 (-3) 


-6.7448 (-4) 




12 


-2.4071 (-3) 


1.8936 (-2) 


1.9108 (-2) 


-4.5379 (-3) 


2.1756 (-4) 




13 


1.0215 (-2) 


1.3014 (-2) 


1.0859 (-2) 


1.9591 (-3) 


6.6418 (-4) 



A=2.0 


\ P 
I \ 
\ 





0.2 


0.5 


0.8 


1.0 







6. 6333 


(-1) 


6.5663 (-1) 


6.2286 (-1) 


5.6597 (-1) 


5.1954 (-1) 




1 


6. 2832 


(-1) 


6.2033 (-1) 


5.8026 (-1) 


5.1359 (-1) 


4.6008 (-1) 




2 


5.6199 


(-1) 


5.5195 (-1) 


5.0190 (-1) 


4.2003 (-1) 


3.5598 (-1) 




3 


4. 7124 


(-1) 


4.5928 (-1) 


3.9995 (-1) 


3.0454 (-1) 


2.3208 (-1) 




4 


3. 6529 


(-1) 


3.5254 (-1) 


2.8937 (-1) 


1.8915 (-1) 


1.1540 (-1) 




5 


2. 5447 


(-1) 


2.4285 (-1) 


1.8488 (-1) 


9.2889 (-2) 


2.7138 (-2) 




6 


1. 4893 


(-1) 


1.4073 (-1) 


9.8168 (-2) 


2.7404 (-2) 


-2.2682 (-2) 




7 


5. 7413 


(-2) 


5.4681 (-2) 


3.6001 (-2) 


-4.9330 (-3) 


-3.6350 (-2) 




8 


-1. 3698 


(-2) 


-9.7870 (-3) 


-3.7792 (-4) 


-1.0030 (-2) 


-2.5404 (-2) 




9 


-6. 1018 


(-2) 


-5.0573 (-2) 


-1.4595 (-2) 


1.1624 (-3) 


-4.9948 (-3) 




10 


-8. 4305 


(-2) 


-6.8873 (-2) 


-1.3479 (-2) 


1.6821 (-2) 


1.1925 (-2) 




11 


-8. 6195 


(-2) 


-6.8555 (-2) 


-4.9232 (-3) 


2.7909 (-2) 


1.8444 (-2) 




12 


-7. 1549 


(-2) 


-5.5188 (-2) 


4.1327 (-3) 


3.0222 (-2) 


1.4479 (-2) 




13 


-4. 6553 


(-2) 


-3.4932 (-2) 


9.1247 (-3) 


2.4418 (-2) 


4.7874 (-3) 




\, 


















1.2 




1.5 


2.0 


5.0 


10.0 







4. 6970 


(-1) 


3.9438 (-1) 


2.8271 (-1) 


4.1460 (-2) 


5.9847 (-3) 




1 


4. 0376 


(-1) 


3.2134 (-1) 


2.0696 (-1) 


1.5787 (-2) 


1.1832 (-3) 




2 


2.9065 


(-1) 


2.0014 (-1) 


8.9034 (-2) 


-1.1646 (-2) 


-2.6407 (-3) 




3 


1.6107 


(-1) 


6.9906 (-2) 


-2.1766 (-2) 


-1.7831 (-2) 


-1.6583 (-3) 




4 


4. 6708 


(-2) 


-3.2114 (-2) 


-8.5503 (-2) 


-3.1794 (-3) 


1.4056 (-3) 




5 


-3.0151 


(-2) 


-8.4575 (-2) 


-8.9692 (-2) 


1.1935 (-2) 


1.8256 (-3) 




6 


-6. 2584 


(-2) 


-8.7743 (-2) 


-5.0978 (-2) 


1.0892 (-2) 


-5.0880 (-4) 




7 


-5.8611 


-2) 


-5.9240 (-2) 


-1.3490 (-3) 


-2.4198 (-3) 


-1.7496 (-3) 




8 


-3. 5329 


(-2) 


-2.3059 (-2) 


3.1301 (-2) 


-1.1062 (-2) 


-2.0353 (-4) 




9 


-1.0687 


-2) 


2.0473 (-3) 


3.6461 (-2) 


-5.7689 (-3) 


1.4767 (-3) 




10 


3. 4017 


-3) 


9.2652 (-3) 


2.1731 (-2) 


5.5909 (-3) 


7.3707 (-4) 




11 


4. 0067 


-3) 


3.2242 (-3) 


3.1477 (-3) 


9.1154 (-3) 


-1.0619 (-3) 




12 


-4. 2982 


-3) 


-5.9353 (-3) 


-6.8639 (-3) 


1.5814 (-3) 


-1.0761 (-3) 




13 


-1.3803 


-2) 


-9.7357 (-3) 


-5.9984 (-3) 


-7.0051 (-3) 


5.6923 (-4) 
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Pi (p, h) 





\, 












ft=5.0 


I \ 





0.2 


0.5 


0.8 


1.0 







1.2202 (-1) 


1.2171 f !> 


1.2033 (-1) 


1.1776 (-1) 


1.1547 (-1) 




1 


L.2083 (-1) 


1.2047 (-1) 


1 L863 (-1) 


1. 1530 (-1) 


1.1235 (-1) 




2 


1. IMS (-1) 


1.1797 (-1) 


1.1528 (-1) 


1.1050 (-1) 


1.0628 (-1) 




3 


1.1502 (-1) 


1. 1427 (-1) 


1. ion (-1) 


1.0358 (-1) 


9.7641 (-2) 




4 


1.1051 (-1) 


L.0947 (-1) 


1.(1117 (-1) 


9.4882 (-2) 


8.6923 (-2) 




5 


L.0503 (-1) 


1.0368 (-1) 


9.6759 (-2) 


8.4810 (-2) 


7.4740 (-2) 




6 


9.8700 (-2) 


9.7005 (-2) 


8 8421 (-2) 


7.3820 (-2) 


6.1761 (-2) 




7 


9.1627 (-2) 


8.9600 (-2) 


7.9414 (-2) 


0.2393 (-2) 


4.8664 (-2) 




8 


8.3940 (-2) 


8.1617 (-2) 


7.0009 (-2) 


5.1001 (-2) 


3.6090 (-2) 




9 


7.5798 (-2) 


7.3219 (-2) 


6.0477 (-2) 


4.0086 (-2) 


2.4596 (-2) 




10 


6.7330 (-2) 


6.4572 (-2) 


5.1080 (-2) 


3.0033 (-2) 


1.4624 (-2) 




11 


5.8690 (-2) 


5.5842 (-2) 


4.2057 (-2) 


2.1148 (-2) 


6.4728 (-3) 




12 


5.0028 (-2) 


4.7193 (-2) 


3.3019 (-2) 


1.3649 (-2) 


2.9255 (-4) 




13 


4.1487 (-2) 


3.8775 (-2) 


2.5939 (-2) 


7.6574 (-3) 


-3.9175 (-3) 




I \ 
\ 


1.2 


1.5 


2.0 


5.0 


10.0 







1.1277 (-1) 


1.0806 (-1) 


9.8893 (-2) 


4.4590 (-2) 


1.1307 (-2) 




1 


1.0889 (-1) 


1.0291 (-1) 


9.1512 (-2) 


3.1724 (-2) 


5.0832 (-3) 




2 


L.0140 (-1) 


9.3079 (-2) 


7.7719 (-2) 


1.1644 (-2) 


-2.2205 (-3) 




3 


9.0832 (-2) 


7.9469 (-2) 


5.9298 (-2) 


-7.14(H) (-3) 


-5.0443 (-3) 




4 


7.7933 (-2) 


6 3284 (-2) 


3.8554 (-2) 


-1.7107 (-2) 


-2.3029 (-3) 




5 


6.3572 (-2) 


L5910 (-2) 


L.7958 (-2) 


-1.6477 (-2) 


2. 1003 (-3) 




6 


L8683 (-2) 


2.8750 (-2) 


-2.3434 (-4) 


-7.1132 (-3) 


3.6874 (-3) 




7 


3.4173 (-2) 


1.3078 (-2) 


-1.4299 (-2) 


4.4049 (-3) 


1.2200 (-3) 




8 


2.0853 (-2) 


-8.7710 (-5) 


-2.3258 (-2) 


1.1755 (-2) 


-2. 1739 (-3) 




9 


9.3700 (-3) 


-1.0080 (-2) 


-2.6958 (-2) 


1.1721 (-2) 


-2.9248 (-3) 




10 


L.6819 (-4) 


-1.6621 (-2) 


-2.6000 (-2) 


5.4042 (-3) 


-5.4759 (-4) 




li 


-6.5366 (-3) 


-1.9815 (-2) 


-2.1555 (-2) 


-2.9324 (-3) 


2. 1682 (-3) 




12 


-1.0753 (-2) 


-2.0093 (-2) 


-1.5101 (-2) 


-8.53(i9 (-3) 


2.3570 (-3) 




13 


-1.2092 (-2) 


-1.8125 (-2) 


-8.1402 (-3) 


-8.7680 (-3) 


4.6573 (-5) 


h=lQ.O 


V 





0.2 


0.5 


0.8 


1.0 







3.1182 (-2) 


3. li ivi (-2) 


3. IOCS (-2) 


3.0890 (-2) 


3.072S (-2) 




1 


3.1105 (-2) 


3.1081 (-2) 


3.0953 (-2) 


3.0718 (-2) 


3.0504 (-2) 
3.0059 (-2) 




2 


3.0951 (-2) 


3.0914 (-2) 


3.0725 (-2) 


:U>376 (-2) 




3 


3.0720 (-2) 


3.0666 (-2) 


3.0384 (-2) 


2.9869 (-2) 


2. '.MOO (-2) 




4 


3.0414 (-2) 


3.0337 (-2) 


2.9936 (-2) 


2.9202 (-2) 


2.8539 (-2) 




5 


3.0035 (-2) 


2.9930 (-2) 


2.9382 (-2) 


2.8385 (-2) 


2.7488 (-2) 




6 


2.9684 (-2) 


2.9446 (-2) 


2.8728 (-2) 


2.7429 (-2) 


2.6266 (-2) 




7 


2.9064 (-2) 


2.8889 (-2) 


2.7981 (-2) 


2.6344 (-2) 


2. L889 (-2) 




8 


2.8470 (-2) 


2.8261 (-2) 


2.7145 (-2) 


2.5145 (-2) 


2.3381 (-2) 




9 


2.7825 (-2) 


2.7566 (-2) 


2.6228 (-2) 


2.3848 (-2) 


2 L764 (-2) 




10 


2.7113 (-2) 


2.6808 (-2) 


2.5239 (-2) 


2.2466 (-2) 


2.0083 (-2) 




11 


2.6343 (-2) 


2.5991 (-2) 


2.4184 (-2) 


2.1018 (-2) 


1.S302 (-2) 




12 


2.5520 (-2) 


2.5119 (-2) 


2.3074 (-2) 


1.9520 (-2) 


L.6506 (-2) 




13 


2.4647 (-2) 


2.4198 (-2) 


2.1916 (-2) 


1.7990 (-2) 


1.4702 (-2) 




\, 














'\ 


1.2 


1.5 


2.0 


5.0 


10.0 







3.0532 (-2) 


3.0176 (-2) 


2.9429 (-2) 


2.2412 (-2) 


1.1118 (-2) 




1 


3.0245 (-2) 


2.9776 (-2) 


2.8801 (-2) 


2.0042 (-2) 


7.8735 (-3) 




2 


2.9077 (-2) 


2.8988 (-2) 


2.7568 (-2) 


1.5688 (-2) 


2.8106 (-3) 




3 


2.8838 (-2) 


2.7832 (-2) 


2.5777 (-2) 


1.0051 (-2) 


-1.9193 (-3) 




4 


2.7746 (-2) 


2.6336 (-2) 


2.3496 (-2) 


4.0196 (-3) 


-4.4729 (-3) 




5 


2.6422 (-2) 


2.4540 (-2) 


2.0810 (-2) 


-1.4972 (-3) 


-4.1608 (-3) 




6 


2.4890 (-2) 


2.2487 (-2) 


1.7818 (-2) 


-5.7256 (-3) 


-1.6836 (-3) 




7 


2.3181 (-2) 


2.0229 (-2) 


1.4627 (-2) 


-8.1599 (-3) 


1.3314 (-3) 




8 


2.1327 (-2) 


1.7821 (-2) 


1.1348 (-2) 


-8.6369 (-3) 


3.2194 (-3) 




9 


1.9360 (-2) 


1.5318 (-2) 


8.0921 (-3) 


-7.3437 (-3) 


3.1141 (-3) 




10 


1.7318 (-2) 


1.2780 (-2) 


4.9635 (-3) 


-4.7585 (-3) 


1.3036 (-3) 




11 


1.5235 (-2) 


1.0261 (-2) 


2.0567 (-3) 


-1.5409 (-3) 


-1.0422 (-3) 




12 


1.3148 (-2) 


7.8155 (-3) 


-5.4776 (-4) 


1.6059 (-3) 


-2.5824 (-3) 




13 


1.1089 (-2) 


5.4916 (-3) 


-2.7864 (-3) 


4.0706 (-3) 


-2.5433 (-3) 



*The figures in parentheses in tables 1, 2, and 3 indicate the power of 10 by which the adjacent entry is to be multiplied; e.g., 2.2831 (—2) =0.022831. 



source -plane, table 1 may be used in two different 
ways: 

(a) If the angular distribution data is already in the 
form of gi coefficients of a Legendre expansion as 
given in (4), a new set of coefficients p'i{p,h) may be 
obtained from the entries in table 1 according 
to 

pl( P) h)={lp l ^(p ) h)+(l+l)p l ^( P} h)} 1(21+1) (32) 

which can be combined with the (ji data as in (6). 



(b) If the source data are given in terms of the 
actual angular distribution g(6), which could be 
either a mathematical expression or experimental 
data, it would be simpler to modify the source 
angular distribution coefficients, that is, to generate 
Legendre coefficients g[ of the quantity g'(d) = 
cos 6 g(6). These gi coefficients could then be com- 
bined directly with p t values from table 1 in the 
same way as has been specified in (6) for an iso- 
tropic detector. 
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PLANE OF S 




K SOURCE, S 



Figure 3. Three examples of detector positions in which the 
distances off-axis and off source-plane, measured in disk 
radii, are (A) p = 0.5, h = 0.5, (B) p=l, h=\ and (C) p = 2, 
h=l. 




2tt 



0.2 



cos 9 



Figure 4. Reconstruction of the response function \p(0,p,h) 
according to (5) for (A) p = .5, h = .5, (B) p=l, h=\ and (C) 
p = 2, h = l using table 1 and summing over terms 0<Z<9. 
The reconstructed functions are shown as dotted lines and the 
exact functions as solid lines. Note that the scale for (B) 
is on the right-hand side. 



2.3. Limitation of the Legendre Expansion Method 

The utility of this method is dependent on the 
oscillatory nature of P z (cos#). Within any finite 
region between cos 6 and cos 0+ Acos 6 the increasing 
number of these oscillations, for increasing Z, yields 
decreasing values of p z and g z as a result of higher 
and higher order cancellations. If both functions 
#s(0) and g(6) are concentrated mostly within very 
narrow regions of cos 0, adequate convergence of (6) 
may not be achieved without going to exorbitantly 
high Z's. 

Evaluation of the unscattered radiation com- 
ponent provides an example of the above limita- 
tion. This problem has been qualitatively discussed 
in section 1, and will be treated in some detail in 
the following section. 



3. Unscattered Radiation Component 

3.1. t>1: Legendre Expansion 

As in the rectangular source case [3], the Legendre 
sum (6) converges adequately for an angular dis- 
tribution of the form (1) when the barrier thickness 
T=\d is of the order of one or greater. In this 
region the behavior of g(6) at the barrier exit surface 
lies between that of a "perfectly diffuse" plane 
source {g(6) = constant} and that of a plane normal 
source {g(d) cc8(6— 0)} . The g° coefficients are 
easily obtained from appropriate combinations of 
tabulated exponential integrals [17], and a sample 
tabulation of gfs is given in reference [3] for 
0<Z<13 over the range 0.001 <r< 100.0. 



3.2. 0<r<l; A>0: Power Series Expansion 

Since one of the chief reasons for separate treat- 
ment of the unscattered radiation component is the 
inadequate convergence of (6) for r<<0> an al- 
ternative means of evaluating the integral 



I°(r }Pj h)^j s (dS/r') exp (-r/cos 6) 



(33) 



is desirable. Again, as in [3], no success was had 
in obtaining a closed -form result for I°(r,p,h) in (33). 
However, Sievert [18] has shown that by expanding 
exp (— r/cos 6) as a (r/cos 6) n power series, a solution 
corresponding to eq (10) in reference [3] may 
be written 



47T w =o 



(34) 



in which the q n (p,h) coefficients are evaluations 
of the integrals 

^^=wf >+ K2)f ~' 2 cos ~ i (W) RdB 

(35) 
for l<p< °° and 



q_ n (p,h)-- 



n\h n \J 



(h 2 +R 2 ) 2 2wBdB 



+J i 1+ >+£') ,t ~ 1 2 cos- 1 ( y+ 2 ^ 1 )^g\ (36) 

for 0<p<l. As in the Pi(p,h) evaluation, three 
distinct expressions arise from (35) and (36) de- 
pending on whether n is zero, even or odd, but by 
proper manipulation the distinctions between re- 
sults from (35) and (36) were removed. In reference 
[18], q n (p,h) coefficients are given only for n=0, 1, 
and 2, of which qi(p,h) is not exact. 

The zeroth coefficient q Q (p,h), when multiplied by 
the source strength factor <r/47r, is the radiation 
field from a disk source in a nonattenuating medium 
(i.e., r = 0). This can be written, for A>0, as [18] 
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^ (p,/0=/ T °=o47r/cr=7rln{(l+A 2 -p 2 

+V(l+/^ 2 -p 2 ) 2 +4p% 2 )/2A 2 }. (37) 

Higher even-n terms (i.e., n(e)=2, 4, 6, ... ) are 
closed-form purely algebraic expressions, the next two 
of which are 



&(p ,&)= 



IT 

= 2!/T 2 



and 



S4(p,A)= i j^(A 2 + P 2 + i). 



(38) 



(39) 



These are included in the general expression for the 
even-rt terms: 



TflT (?)(©( 



A 2 +P 2 +l 



-Y 
-1/ 



p 2 +i 



,.«<M4-l/2) 

-£(-i) 

m=0 V "* 



The odd-n terms (i.e., «(o) = l, 3, 5, 
included in the series-expansion solution 



A 2 +P 2 - 
(2i\(_p V'T/ __^ \" /2 

M/v+i/ Lw+p*+i/ 

. ) are all 



q««»(p,h)-- 



+ 



A 2 + P 2 



p 2 + 



T{sG0OD( 

l&wv+v LU 2 +p 2 +i/ 



a 2 +p 2 +a-. 



2! 

■S(- 

ra=0 



Cf)(p^)1}}. <«> 



which strongly resembles (40), except that here the 
quantities (%?) and ( w 4 2 ) are now coefficients of the 
infinite series expansion of half-integral powers of 
a binomial [19], i.e., 

(l+xy»=±( n P)x m . (42) 

These are readily generated using the recursion 

CSV 

(n\T 



-\ 



o 



=n/2 ) and 

2(m— 1)— n 
2 /// 



/ n/2 \ 

U-i/J 



(43) 



which is also valid when 7i/2 is an integer, since ( n m), 
using (43), is zero for m^>nl2. 

By setting p=0 in (37), (40) and (41) and substi- 
tuting the resulting values for q n (pjh) in (34), we 
find the on-axis unscattered flux to be 



P(0,h) = (cr/4T)27r\ In 



[- 



Va 2 +i 



+s 



(-1) 



" {(^y-'H <"> 



1 71 • 7Z! 

This series can be identified with the familiar on- 
axis solution in terms of exponential integrals [20] 

I°(Q,h) = (a/4T) I V + ^- 2*rdr 

= (<r/4ir)2ir{S 1 (r) - E x (r^'hJ+i/h) } . (45) 

Equations (37), (40) and (41) were programed 
for evaluation on an IBM 704 computer and a 
resulting set of q n (p,h) coefficients for 0<p<10, 
0.1 </i<10 and 0<n<9 is given in table 2. 



3.3. Comparison of Legendre and Power Series 
Expansions 

As a partial check on tables 1 and 2, and because 
the unscattered radiation component is of some gen- 
eral interest, the quantities / (t,p,A)47t/o" as defined 
in (33) were evaluated using the methods of the 
preceding two sections, and are given in table 3. 
Numbers preceded by asterisks were evaluated using 
the Legendre expansion approach (6) with g\ values 
from reference [3] and 0<Z<13; the others are 
evaluated using the r n series (34) for 0<n<20. 

Some notion of the convergence for each method 
may be obtained from table 4, in which partial sums 
over 0</<13 and 0<^<13 for r = 1.0 are given for 
the three situations shown in figure 3. A qualitative 
comment about the Legendre sums for these three 
cases is that in each case the first two terms con- 
tribute over 80 percent of the result. 

For scattered radiation the tendency is toward 
even faster convergence. For example, the Fermi 
approximation [21] 



0(0) = (l + V3cos0)/(l+V'3/2), 



(46) 



used in reference [7] for the angular distribution of 
neutrons from a reactor thermal-column face, has 
been shown by Placzek [22] to be in error by less 
than one percent for 0.1<cos 0<O> with a maximum 
error of 7.2 percent at cos 0=0. The Legendre ex- 
pansion coefficients for the g(6) in (46) are: ^0= 
1.072, #! = 0.6188 and gi>2=0, which may be used 
with the pi values in table 1 according to (6). 
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Table 2. Coefficients q r n (p,h) [ = h n q n (p,h)] of the (r//t) n power series expansion (34) of the unscattered component Uirlo)I°(T,p,h) 
of the radiation flux, from a circular plane isotropic source, incident on a detector at position p,h separated from the source- 
plane by an attenuating layer of thickness r=ixt. Each term contains the factor h n in order to include meaningful tabulations 
in the region p^>l,h^O as in {47). 

q n ' (p,h) =h n g n (p,h) 



h=0 


X 
X p 
n \ 
X 





0.2 


0.5 


0.8 


1.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 




[No physical solution within the radiating surface.] 






n \ 
\ 


1.2 


1.5 


2.0 


5.0 


10.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


3.6157 (0) 
-2.9533 (0) 

1.5708 (0) 
-6.8470 (-1) 

2.5395 (-1) 
-8.1621 (-2) 

2.3069 (-2) 
-5.8064 (-3) 

1.3158 (-3) 
-2.7095 (-4) 


1.8465 (0) 
-2.2364 (0) 

1.5708 (0) 
-8.2992 (-1) 

3.5997 (-1) 
-1.3336 (-1) 

4.3179 (-2) 
-1.2414 (-2) 

3.2080 (-3) 
-7.5263 (-4) 


9.0378 (-1) 
-1.6252 (0) 

1.5708 (0) 
-1.0803 (0) 

5.8905 (-1) 
-2.6896 (-1) 

1.0617 (-1) 
-3.6999 (-2) 

1.1551 (-2) 
-3.2677 (-3) 


1.2825 (-1) 
-6.3151 (-1) 

1.5708 (0) 
-2.6311 (0) 

3.3379 (0) 
-3.4200 (0) 

2.9467 (0) 
-2.1951 (0) 

1.4425 (0) 
-8.4902 (-1) 


3.1574 (-2) 
-3.1455 (-1) 

1.5708 (0) 
-5.2425 (0) 

1.3155 (1) 
-2.6475 (1) 

4.4507 (1) 
-6.4288 (1) 

8.1446 (1) 
-9.1932 (1) 


h=0.l 


\ 

\ P 
n \ 
\ 





0.2 


0.5 


0.8 


1.0 





1 
2 
3 

4 
5 
6 
7 
8 
9 


1.4499 (1) 
-5.6862 (0) 

1.5708 (0) 
-3.5397 (-1) 

6.6759 (-2) 
-1.0736 (-2) 

1.4985 (-3) 
-1.8441 (-4) 

2.0270 (-5) 
-2.0120 (-6 


1.4373 (1) 
-5.6239 (0) 

1.5708 (0) 
-3.6436 (-1) 

7.1995 (-2) 
-1.2326 (-2) 

1.8581 (-3) 
-2.5003 (-4) 

3.0378 (-5) 
-3.3635 (-6) 


1.3619 (1) 
-5.2804 (0) 

1.5708 (0) 
-4.1803 (-1) 

9.9484 (-2) 
-2.1058 (-2) 

3.9747 (-3) 
-6.7332 (-4) 

1.0316 (-4) 
-1.4403 (-5) 


1.1476 (1) 
-4.5455 (0) 

1.5708 (0) 
-5.1342 (-1) 

1.5054 (-1) 
-3.8949 (-2) 

8.9266 (-3) 
-1.8285 (-3) 

3.3788 (-4) 
-5.6805 (-5) 


7.3908 (0) 
-3.7113 (0) 

1.5708 (0) 
-5.9576 (-1) 

1.9766 (-1) 
-5.7335 (-2) 

1.4676 (-2) 
-3.3511 (-3) 

6.8960 (-4) 
-1.2903 (-4) 




\ 
\ P 
n X 

X 


1.2 


1.5 


2.0 


5.0 


10.0 






1 

2 
3 

4 
5 
6 
7 
8 
9 


3.6033 (0) 
-2.9136 (0) 

1.5708 (0) 
-6.8715 (-1) 

2.5525 (-1) 
-8.2135 (-2) 

2.3238 (-2) 
-5.8551 (-3) 

1.3282 (-3) 
-2.7379 (-4) 


1.8269 (0) 
-2.2276 (0) 

1.5708 (0) 
-8.3178 (-1) 

3.6128 (-1) 
-1.3399 (-1) 

4.3419 (-2) 
-1.2493 (-2) 

3.2312 (-3) 
-7.5869 (-4) 


9.0031 (-1) 
-1.6225 (0) 

1.5708 (0) 
-1.0816 (0) 

5.9036 (-1) 
-2.6977 (-1) 

1.0657 (-1) 
-3.7159 (-2) 

1.1608 (-2) 
-3.2857 (-3) 


1.2819 (-1) 
-6.3138 (-1) 

1.5708 (0) 
-2.6316 (0) 

3.3393 (0) 
-3.4220 (0) 

2.9489 (0) 
-2.1971 (0) 

1.4440 (0) 
-8.5009 (-1) 


3.1571 (-2) 
-3.1454 (-1) 

1.5708 (0) 
-5.2428 (0) 

1.3157 (1) 
-2.6479 (1) 

4.4516 (1) 
-6.4304 (1) 

8.1470 (1) 
-9.1963 (1) 


h=0.2 


X 

n \ 
X 





0.2 


0.5 


0.8 


1.0 






1 

2 
3 
4 
5 
6 
7 
8 
9 


1.0236 (1) 
-5.1510 (0) 

1.5708 (0) 
-3.6743 (-1) 

7.0686 (-2) 
-1.1547 (-2) 

1.6360 (-3) 
-2.0430 (-4) 

2.2788 (-5) 
-2.2952 (-6) 


1.0117 (1) 
-5.0914 (0) 

1.5708 (0) 
-3.7767 (-1) 

7.5922 (-2) 
-1.3161 (-2) 

2.0060 (-3) 
-2.7278 (-4) 

3.3482 (-5) 
-3.7445 (-6) 


9.4210 (0) 
-4.7664 (0) 

1.5708 (0) 
-4.3050 (-1) 

1.0341 (-1) 
-2.2013 (-2 

4.1776 (-3) 
-7.1178 (-4) 

1.0971 (-4) 
-1.5412 (-5) 


7.6461 (0) 
-4.1033 (0) 

1.5708 (0) 
-5.2415 (-1) 

1.5446 (-1) 
-4.0116 (-2) 

9.2316 (-3) 
-1.8991 (-3) 

3.5247 (-4) 
-5.9523 (-5) 


5.3698 (0) 
-3.4555 (0) 

1.5708 (0) 
-6.0468 (-1) 

2.0159 (-1) 
-5.8685 (-2) 

1.5075 (-2) 
-3.4547 (-3) 

7.1350 (-4) 
-1.3399 (-4) 




X 
X p 
n X 


1.2 


1.5 


2.0 


5.0 


10.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


3.2477 (0) 
-2.8313 (0) 

1.5708 (0) 
-6.9432 (-1) 

2.5918 (-1) 
-8.3689 (-2) 

2.3753 (-2) 
-6.0032 (-3) 

1.3660 (-3) 
-2.8243 (-4) 


1.7714 (0) 
-2.2021 (0) 

1.5708 (0) 
-8.3732 (-1) 

3.6521 (-1) 
-1.3586 (-1) 

4.4146 (-2) 
-1.2734 (-2) 

3.3015 (-3) 
-7.7708 (-4) 


8.9009 (-1) 
-1.6145 (0) 

1.5708 (0) 
-1.0857 (0) 

5.9429 (-1) 
-2.7221 (-1) 

1.0775 (-1) 
-3.7643 (-2) 

1.1780 (-2) 
-3.3402 (-3) 


1.2803 (-1) 
-6.3098 (-1) 

1.5708 (0) 
-2.6332 (0) 

3.3432 (0) 
-3.4279 (0) 

2.9556 (0) 
-2.2032 (0) 

1.4488 (0) 
-8.5330 (-1) 


3.1561 (-2) 
-3.1449 (-1) 

1.5708 (0) 
-5.2436 (0) 

1.3161 (1) 
-2.6490 (1) 

4.4542 (1) 
-6.4351 (1) 

8.1542 (1) 
-9.2057 (1) 



258 



Table 2. Coefficients q f n (p,h) [ = h n q n (p,h)] of the (r//i) w power series expansion (34) of the unscattered component (4Tr/(r)I (T t p,h) 
of the radiation flux, from a circular plane isotropic source, incident on a detector at position p,h separated from the source- 
plane by an attenuating layer of thickness r = nt. Each term contains the factor h n in order to include meaningful tabulations 
in the region py>l,h~0 as in (47) — Continued. 

q n '(p,h)=h"q n (p,h) 



h=0.5 







0.2 


0.5 


0.8 


1.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


5.0562 (0) 
-3.8832 (0) 

1.5708 (0) 
-4.4420 (-1) 

9.8175 (-2) 
-1.7967 (-2) 

2.8180 (-3) 
-3.8751 (-4) 

4.7480 (-5) 
-5.2476 (-6) 


L9763 (0) 
-3.83,s;-i (ID 

1.5708 (0) 
-4.5354 (-1) 

1.0341 (-1) 
-1.9734 (-2) 

3.2613 (-3) 
-4.7714 (-4) 

6.2749 (-5) 
-7.5023 (-6) 


4.5353 (0) 
-3.6032 (0) 

1.5708 (0) 
-5.0155 (-1) 

1.3090 (-1) 
-2.9378 (-2) 

5.8178 (-3) 
-1.0318 (-3) 

1.6557 (-4) 
-2.4234 (-5) 


3.6855 (0) 
-3.1826 (0) 

1.57C8 (0) 
-5.8610 (-1) 

1.8195 (-1) 
-4.8879 (-2) 

1.1587 (-2) 
-2.4543 (-3) 

4.6921 (-4) 
-8.1650 (-5) 


2.9550 (0) 
-2.8363 (0) 

1.5708 (0) 
-6.5868 (-1) 

2.2907 (-1) 
-6.6468 (-2) 

1.8090 (-2) 
-4.2497 (-3) 

8.9969 (-4) 
-1.7321 (-4) 




\ 

\ P 
n \ 
\ 


1.2 


1.5 


2.0 


5.0 


10.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


2.2550 (0) 
-2.4916 (0) 

1.5708 (0) 
-7.4054 (-1) 

2.8667 (-1) 
-9.4997 (-2) 

2.7574 (-2) 
-7.1192 (-3) 

1.6543 (-3) 
-3.4929 (-4) 


1.4874 (0) 
-2.0556 (0) 

1.5708 (0) 
-8.7450 (-1) 

3.9270 (-1) 
-1.4935 (-1) 

4.9451 (-2) 
-1.4516 (-2) 

3.8277 (-3) 
-9.1603 (-4) 


8.2602 (-1) 
-1.5625 (0) 

1.5708 (0) 
-1.1135 (0) 

6.2177 (-1) 
-2.8953 (-1) 

1.1626 (-1) 
-4.1153 (-2) 

1.3040 (-2) 
-3.7422 (-3) 


1.2690 (-1) 
-6.2825 (-1) 

1.5708 (0) 
-2.6442 (0) 

3.3707 (0) 
-3.4695 (0) 

3.0026 (0) 
-2.2463 (0) 

1.4823 (0) 
-8.7601 (-1) 


3.1494 (-2) 
-3.1416 (-1) 

1.5708 (0) 
-5.2491 (0) 

1.3188 (1) 
-2.6573 (1) 

4.4727 (1) 
-6.4683 (1) 

8.2044 (1) 
-9.2716 (1) 


h=l.Q 


\ 
\ P 
n \ 
\ 


o 


0.2 


0.5 


0.8 


1.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


2.1776 (0) 
-2.6026 (0) 

1.5708 (0) 
-6.3824 (-1) 

1.9635 (-1) 
-4.8767 (-2) 

1.0181 (-2) 
-1.8368 (-3) 

2.9219 (-4) 
-4.1608 (-5) 


2.1463 (0) 
-2.5805 (0) 

1.5708 (0) 
-6.4562 (-1) 

2.0159 (-1) 
-5.0998 (-2) 

1.0886 (-2) 
-2.0161 (-3) 

3.3053 (-4) 
-4.8702 (-5) 


1.9878 (0) 
-2.4687 (0) 

1.5708 (0) 
-6.8363 (-1) 

2.2907 (-1) 
-6.3027 (-2) 

1.4817 (-2) 
-3.0566 (-3) 

5.6368 (-4) 
-9.4190 (-5) 


1.7219 (0) 
-2.2805 (0) 

1.5708 (0) 
-7.5100 (-1) 

2.8013 (0) 
-8.6736 (-2) 

2.3138 (-2) 
-5.4471 (-3) 

1.1504 (-3) 
-2.2055 (-4) 


1.5118 (0) 
-2.1288 (0) 

1.5708 (0) 
-8.0972 (-1) 

3.2725 (-1) 
-1.1015 (-1) 

3.1998 (-2) 
-8.2091 (-3) 

1.8895 (-3) 
-3.9469 (-4) 




\ 

n \ 
\ 


1.2 


1.5 


2.0 


5.0 


10.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


1.2993 (0) 
-1.9698 (0) 

1.5708 (0) 
-8.7743 (-1) 

3.8485 (-1) 
-1.4065 (-1) 

4.4362 (-2) 
-1.2351 (-2) 

3.0830 (-3) 
-6.9802 (-4) 


1.0123 (0) 

-1.7394 (0) 

1.5708 (0) 
-9.9202 (-1) 

4.9087 (-1) 
-2.0193 (-1) 

7.1540 (-2) 
-2.2333 (-2) 

6.2418 (-3) 
-1.5806 (-3) 


6.6582 (-1) 
-1.4173 (0) 

1.5708 (0) 
-1.2063 (0) 

7.1995 (-1) 
-3.5482 (-1) 

1.4981 (-1) 
-5.5516 (-2) 

1.8369 (-2) 
-5.4956 (-3) 


1.2302 (-1) 
-6.1877 (-1) 

1.5708 (0) 
-2.6832 (0) 

3.4688 (0) 
-3.6193 (0) 

3.1736 (0) 
-2.4045 (0) 

1.6064 (0) 
-9.6077 (-1) 


3.1257 (-2) 
-3.1298 (-1) 

1.5708 (0) 
-5.2687 (0) 

1.3286 (1) 
-2.6869 (1) 

4.5389 (1) 
-6.5876 (1) 

8.3854 (1) 
-9.5095 (1) 


h=2.0 


V 

n \ 
\ 





0.2 


0.5 


0.8 


1.0 





1 
2 
3 
4 
5 
6 
7 
8 
9 


7.0103 (-1) 
-1.4833 (0) 

1.5708 (0) 
-1.1102 (0) 

5.8905 (-1) 
-2.5030 (-1) 

8.8721 (-2) 
-2.6983 (-2) 

7.1878 (-3) 
-1.7037 (-3) 


6.9603 (-1) 
-1.4777 (0) 

1.5708 (0) 
-1.1148 (0) 

5.9428 (-1) 
-2.5382 (-1) 

9.0473 (-2) 
-2.7684 (-2) 

7.4236 (-3) 
-1.7722 (-3) 


6.7068 (-1) 
-1.4491 (0) 

1.5708 (0) 
-1.1391 (0) 

6.2177 (-1) 
-2.7251 (-1) 

9.9902 (-2) 
-3.1511 (-2) 

8.7303 (-3) 
-2.1583 (-3) 


6.2738 (-1) 
-1.3995 (0) 

1.5708 (0) 
-1.1831 (0) 

6.7282 (-1) 
-3.0818 (-1) 

1.1843 (-1) 
-3.9278 (-2) 

1.1475 (-2) 
-2.9996 (-3) 


5.9133 (-1) 
-1.3572 (0) 

1.5708 (0) 
-1.2226 (0) 

7.1995 (-1) 
-3.4218 (-1) 

1.3672 (-1) 
-4.7227 (-2) 

1.4395 (-2) 
-3.9319 (-3) 




\ 
\ P 
n \ 
\ 


1.2 


1.5 


2.0 


5.0 


10.0 





1 

2 
3 
4 
5 
6 
7 
8 
9 


5.5179 (-1) 
-1.3098 (0) 

1.5708 (0) 
-1.2693 (0) 

7.7754 (-1) 
-3.8510 (-1) 

1.6060 (-1) 
-5.7991 (-2) 

1.8501 (-2) 
-5.2954 (-3) 


4.8995 (-1) 
-1.2327 (0) 

1.5708 (0) 
-1.3517 (0) 

8.8357 (-1) 
-4.6789 (-1) 

2.0899 (-1) 
-8.0941 (-2) 

2.7732 (-2) 
-8.5329 (-3) 


3.9168 (-1) 
-1.1012 (0) 

1.5708 (0) 
-1.5157 (0) 

1.1127 (0) 
-6.6249 (-1) 

3.3307 (-1) 
-1.4532 (-1) 

5.6119 (-2) 
-1.9471 (-2) 


1.0970 (-1) 
-5.8485 (-1) 

1.5708 (0) 
-2.8335 (0) 

3.8615 (0) 
-4.2401 (0) 

3.9066 (0) 
-3.1058 (0) 

2.1742 (0) 
-1.3612 (0) 


3.0342 (-2) 
-3.0839 (-1) 

1.5708 (0) 
-5.3464 (0) 

1.3679 (1) 
-2.8063 (1) 

4.8085 (1) 
-7.0780 (1) 

9.1364 (1) 
-1.0506 (2) 
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Table 2. Coefficients q' n (p,h) [ = h n q n (p,h)] of the (r/h) n power series expansion (34) of the unscattered component (4* I a) 1° (r,p y h) 
of the radiation flax, from a circular plane isotropic source, incident on a detector at position p,h separated from the source- 
plane by an attenuating layer of thickness r = p.t. Each term contains the factor h n in order to include meaningful tabulations 
in the regular p^>l,h^O as in (47) — Continued 

q n '(p,h)=h n q n (p,h) 



A = 5.0 


n \ 
\ 





0.2 


0.5 


0.8 


1.0 







1.2322 (-1) 


1.2303 (-1) 


1.2206 (-1) 


1.2031 (-1) 


1.1873 (-1) 




1 


-6.2216 (-1) 


-6.2168 (-1) 


-6.1922 (-1) 


-6.1471 (-1) 


-6.1063 (-1) 




2 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 




3 


-2.6440 (0) 


-2.6461 (0) 


-2.6568 (0) 


-2.6767 (0) 


-2.6949 (0) .J 




4 


3.3379 (0) 


3.3432 (0) 


3.3707 (0) 


3.4217 (0) 


3.4688 (0) 




5 


-3.3713 (0) 


-3.3793 (0) 


-3.4215 (0) 


-3.5003 (0) 


-3. 5735 (0) J 




6 


2.8376 (0) 


2.8467 (0) 


2.8946 (0) 


2.9846 (0) 


3.0689 (0) i 




7 


-2.0473 (0) 


-2.0555 (0) 


-2.0993 (0) 


-2.1820 (0) 


-2.2599 (0) 




8 


1.2925 (0) 


1.2988 (0) 


1.3324 (0) 


1.3962 (0) 


1.4568 (0) 




9 


-7.2531 (-1) 


-7.2949 (-1) 


-7.5174 (-1) 


-7.9429 (-1) 


-8.3502 (-1) 




\ 
\ P 


1.2 


1.5 


2.0 


5.0 


10.0 




n \ 
\ 

















1.1686 (-1) 


1.1355 (-1) 


1.0699 (-1) 


6.2828 (-2) 


2.5193 (-2) 




1 


-6.0574 (-1) 


-5.9704 (-1) 


-5.7943 (-1) 


-4.4372 (-1) 


-2.8111 (-1) 




2 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 




3 


-2.7170 (0) 


-2.7572 (0) 


-2.8422 (0) 


-3.7163 (0) 


-5.8610 (0) 




4 


3.5264 (0) 


3.6325 (0) 


3.8615 (0) 


6.6104 (0) 


1.6428 (1) 




5 


-3.6636 (0) 


-3.8314 (0) 


-4.2021 (0) 


-9.4298 (0) 


-3.6895 (1) 




6 


3.1734 (0) 


3.3702 (0) 


3.8150 (0) 


1.1237 (1) 


6.9160 (1) 




7 


-2.3574 (0) 


-2.5430 (0) 


-2.9723 (0) 


-1.1505 (1) 


-1.1130 (2) 




8 


1.5331 (0) 


1.6802 (0) 


2.0286 (0) 


1.0332 (1) 


1.5695 (2) 




9 


-8.8671 (-1) 


-9.8760 (-1) 


-1.2322 (0) 


-8.2664 (0) 


-1.9705 (2) 


h = 10.0 


\ 
\ P 





0.2 


0.5 


0.8 


1.0 







3.1260 (-2) 


3.1248 (-2) 


3.1183 (-2) 


3.1064 (-2) 


3.0955 (-2) 




1 


-3.1338 (-1) 


-3.1332 (-1) 


-3.1299 (-1) 


-3.1239 (-1) 


-3.1184 (-1) 




2 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 




3 


-5.2491 (0) 


-5.2501 (0) 


-5.2556 (0) 


-5.2657 (0) 


-5.2750 (0) 




4 


1.3155 (1) 


1.3161 (1) 


1.3188 (1) 


1.3239 (1) 


1.3286 (1) 




5 


-2.6377 (1) 


-2.6392 (1) 


-2.6475 (1) 


-2.6630 (1) 


-2.6772 (1) 




6 


4.4071 (1) 


4.4106 (1) 


4.4292 (1) 


4.4637 (1) 


4.4957 (1) 




7 


-6.3116 (1) 


-6.3180 (1) 


-6.3512 (1) 


-6.4133 (1) 


-6.4710 (1) 




8 


7.9093 (1) 


7.9188 (1) 


7.9691 (1) 


8.0629 (1) 


8.1501 (1) 




9 


-8.8102 (1) 


-8.8226 (1) 


-8.8881 (1) 


-9.0106 (1) 


-9.1247 (1) 




\ 

\ P 


1.2 


1.5 


2.0 


5.0 


10.0 




n \ 
\ 

















3.0823 (-2) 


3.0582 (-2) 


3.0074 (-2) 


2.5072 (-2) 


1.5708 (-2) 




1 


-3.1117 (-1) 


-3.0995 (-1) 


-3.0736 (-1) 


-2.8060 (-1) 


-2.2207 (-1) 




2 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 


1.5708 (0) 




3 


-5.2864 (0) 


-5.3073 (0) 


-5.3522 (0) 


-5.8645 (0) 


-7.4117 (0) 




4 


1.3344 (1) 


1.3450 (1) 


1.3679 (1) 


1.6428 (1) 


2.6245 (1) 




5 


-2.6947 (1) 


-2.7270 (1) 


-2.7971 (1) 


-3.6829 (1) 


-7.4395 (1) 




6 


4.5349 (1) 


4.6076 (1) 


4.7666 (1) 


6.8833 (1) 


1.7584 (2) 




7 


-6.5419 (1) 


-6.6735 (1) 


-6.9633 (1) 


-1.1031 (2) 


-3.5647 (2) 




8 


8.2576 (1) 


8.4579 (1) 


8.9016 (1) 


1.5475 (2) 


6.3271 (2) 




9 


-9.2656 (1) 


-9.5289 (1) 


-1.0116 (2) 


-1.9305 (2) 


-9.9885 (2) 



Table 3. Unscattered component (4^\<j)I° (T,p,h) of the radiation flux calculated (a) using table 1 of this work and table 2 of 
reference [3] (numbers preceded by asterisks) and (b) using table 2 of this work (no asterii ' 

(AtcI*) 1° (r, p, h) 



h=0.1 



V 

\ 

\ 





0.01 
0.02 
0.05 

0.1 



0.5 
1.0 
2.0 
5.0 
10.0 



1.45 (1) 

1.39 (1) 

1.34 (1) 

1.20 (1) 

1.01 (1) 

*1.00 (1) 



7.38 
*7.27 



*3.50 (0) 
*1.38 (0) 
*3.08 (-1) 
*7.22 (-3) 
*2.61 (-5) 



1.44 (1) 

1.38 (1) 

1.33 (1) 
1.19 (1) 

1.00 (1) 
*9.98 (0) 

7.34 (0) 
*7.24 (0) 

*3.50 (0) 

*1.38 (0) 

*3.08 (-1) 

*7.22 (-3) 

*2.61 (-5) 



1.36 (1) 

1.31 (1) 

1.26 (1) 

1.13 (1) 

9.57 (0) 
*9.56 (0) 

7.10 (0) 

*7.04 (0) 

*3.46 (0) 
*1.38 (0) 
*3.08 (-1) 
*7.21 (-3) 
*2.61 (-5) 



1.15 (1) 

1.10 (1) 
1.06 (1) 
9.54 (0) 

8.11 (0) 

*8.02 (0) 

6.12 (0) 
*6.07 (0) 

*3.17 (0) 
*1.32 (0) 
*3.05 (-1) 
*7.21 (-3) 
*2.61 (-5) 



7.39 (0) 
7.03 (0) 
6.71 (0) 
5.86 (0) 

4.81 (0) 
*4.72 (0) 

3.45 (0) 
*3.41 (0) 

♦1.66 (0) 
*6.63 (-1) 
*1.49 (-1) 
*3.54 (-3) 
*1.29 (-5) 



3.60 (0) 

3.33 (0) 

3.08 (0) 

2.47 (0) 

1.76 (0) 

*1.65 (0) 

9.94 (-1) 

*9.38 (-1) 

*2.51 (-1) 

*4.25 (-2) 

*2. 11 (-3) 



1.83 (0) 

1.62 (0) 

1.44 (0) 

1.02 (0) 

5.99 (-1) 



*2.09 (-1) 

*2.74 (-2) 
*1.03 (-3) 



9.00 (-1) 
7.53 (-1) 
6.31 (-1) 
3.76 (-1) 

1.66 (-1) 



*4.56 (-2) 

*3.45 (-3) 
*2.49 (-4) 



1.28 (-1) 
7.84 (-2) 
4.81 (-2) 
1.13 (-2) 



*1.36 (-3) 
*2. 17 (-5) 



3.16 (-#) 

1.17 (-2) 
4.31 (-3) 



*1.34 (-4) 
•1.72 (-5) 
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Table 3. Unscattered component (4w/<r)I°(T,p t h) of the radiation flux calculated (a) using table 1 of this work and table 2 of 
reference [8] (numbers preceded by asterisks) and (b) using table 2 of this work (no asterisks) — Continued 

(4W«r) 7° (r, p, h) 



h = 0.2 


\ 

\ 

\ 
\ 





0.2 


0.5 


0.8 


1.0 


1.2 


1.6 


2.0 


5.0 


10.0 






0.01 

0.02 

0.05 

0.1 

0.2 
0.5 

1.0 

2.0 
5.0 
10.0 


1.02 (1) 
9.98 (0) 
9.74 (0) 
9.04 (0) 
8.01 (0) 

6.35 (0) 
*6.39 (0) 

3.37 (0) 
*3.36 (0) 

*1.37 (0) 
*3.08 (-1) 
*7.22 (-3) 
*2.61 (-5) 


1.01 (1) 
9.87 (0) 
9.62 (0) 
8.94 (0) 
7.92 (0) 

6.28 (0) 
*6.36 (0) 

3.35 (0) 
*3.34 (0) 

*1.37 (0) 
*3.08 (-1) 
*7.22 (-3) 
*2.61 (-5) 


9.42 (0) 
9.19 (0) 
8.96 (0) 
8.32 (0) 
7.38 (0) 

5.88 (0) 
*5.88 (0) 

3.19 (0) 
*3.19 (0) 

*1.34 (0) 
*3.06 (-1) 
*7.21 (-3) 
*2.61 (-5) 


7.65 (0) 
7.44 (0) 
7.25 (0) 
6.71 (0) 
5.93 (0) 

4.71 (0) 
*4.70 (0) 

2.60 (0) 
*2.58 (0) 

*1. 14 (0) 

*2.80 (-1) 
*7.08 (-3) 
*2.61 (-5) 


5.37 (0) 
5.20 (0) 
5.04 (0) 
4.60 (0) 
3.97 (0) 

3.17 (0) 
*3.02 (0) 

1.55 (0) 
*1.55 (0) 

*6.35 (-1) 
*1.42 (-1) 
*3.48 (-3) 
*1.27 (-5) 


3.25 (0) 
3.11 (0) 
2.98 (0) 
2.63 (0) 
2.15 (0) 

1.49 (0) 
•1.46 (0) 


1.77 (0) 
1.67 (0) 
1.57 (0) 
1.31 (0) 

9.78 (-1) 

5.66 (-1) 
*5.36 (-1) 


8.90 (-1) 
8.13 (-1) 
7.43 (-1) 
5.70 (-1) 
3.70 (-1) 

1.62 (-1) 
*1.38 (-1) 


1.28 (-1) 
1.00 (-1) 
7.83 (-2) 
3.76 (-2) 
1.12 (-2) 


3. 16 (-2) 
1.92 (-2) 
1. 17 (-2) 
2.63 (-3) 
1.01 (-3) 




*3.48 (-3) 


*3.16 (-4) 




*5.79 (-1) 

•1.67 (-1) 
*2.06 (-2) 
*1.04 (-4) 
*2.71 (-8) 


"1.37 (-1) 

•1.73 (-2) 
*5.93 (-4) 
*7.67 (-7) 


*2.19 (-2) 
*1.02 (-3) 


*1.62 (-4) 


*4.90 (-6) 


































h=0.5 


\ 

\ 

\ 

\ 





0.2 


0.5 


0.8 


1.0 


1.2 


15 


2.0 


5.0 


10.0 






0.01 

0.02 

0.05 

0.1 

0.2 

0.5 
1.0 

2.0 

5.0 

Ki.o 


5.06 (0) 
4.98 (0) 
4.90 (0) 
4.68 (0) 
4.34 (0) 
3.73 (0) 

2.38 (0) 
*2.38 (0) 

1.16 (0) 
*1. 16 (0) 

*2.93 (-1) 
*7.21 (-3) 
*2.61 (-5) 


4.98 (0) 
4.90 (0) 
4.82 (0) 
4.61 (0) 
4.27 (0) 
3.66 (0) 

2.34 (0) 
*2.34 (0) 

1.14 (0) 
*1. 14 (0) 

*2.90 (-1) 
*7.20 (-3) 
*2.61 (-5) 


4.53 (0) 
4.46 (0) 
4.39 (0) 
4.19 (0) 
3.87 (0) 
3.32 (0) 

2.11 (0) 
*2.11 (0) 

1.03 (0) 
*1.03 (0) 

*2.09 (-1) 
*7.04 (-3) 
*2.61 (-5) 


3.60 (0) 

3.62 (0) 
3.56 (0) 
3.38 (0) 
3.11 (0) 

2.63 (0) 

1.63 (0) 
*1.63 (0) 

7.78 (-1) 
*7.77 (-1) 

*2.03 (-1) 
*5.76 (-3) 
*2.31 (-5) 


2.96 (0) 
2.90 (0) 
2.84 (0) 
2.69 (0) 
2.45 (0) 
2.03 (0) 

1.21 (0) 
•1.21 (0) 

5.43 (-1) 
*5.43 (-1) 

*1.32 (-1) 
*3.29 (-3) 
•1.22 (-5) 


2.26 (0) 
2.21 (0) 
2.16 (0) 
2.02 (0) 
1.81 (0) 
1.47 (0) 

8.07 (-1) 
*8.07 (-1) 

3.26 (-1) 
*3.26 (-1) 

♦6.64 (-2) 
*1.09 (-3) 
*2.29 (-6) 


1.49 (0) 
1.45 (0) 
1.41 (0) 
L.30 (0) 
1.13 (0) 
8.69 (-1) 

4.09 (-1) 
*4.08 (-1) 


8.26 (-1) 
7.95 (-1) 
7.66 (-1) 
6.84 (-1) 
5.68 (-1) 
3.94 (-1) 

1.38 (-1) 
•1.36 (-1) 


1.27 (-1) 
1.15 (-1) 
1.04 (-1) 
7.74 (-2) 
4.74 (-2) 
1.79 (-2) 


3.15 (-2) 
2.58 (-2) 
2.11 (-2) 
1 Hi (-2) 
4.30 (-3) 








*1.72 (-3) 


6.69 (-5) 




*1.30 (-1) 

*1.73 (-2) 
*9.00 (-5) 
*2.52 (-8) 


*2.73 (-2) 

•1.62 (-3) 
*6.03 (-7) 


*1. 87 (-5) 




























h = l.O 


\ 
\ P 

\ 

\ 

\ 





0.2 


0.5 


0.8 


1.0 


1.2 


1.5 


2.0 


5.0 


10.0 






0.01 
0.02 
0.05 

(i. 1 
0.2 

0.5 

1.0 

2.0 

5.0 
10.0 


2.18 (0) 
2.15 (0) 
2.13 (0) 
2.05 (0) 
1.93 (0) 
1.71 (0) 

1.20 (0) 
♦1.21 (0) 

6.64 (-1) 
*6.64 (-1) 

2.05 (-1) 
*2.06 (-1) 

*6.54 (-3) 
*2.58 (-5) 


2.15 (0) 
2.12 (0) 
2.10 (0) 
2.02 (0) 
1.90 (0) 
1.69 (0) 

1.18 (0) 
*1. 18 (0) 

6.51 (-1) 
*6.51 (-1) 

2.01 (-1) 
*2.01 (-1) 

*6.42 (-3) 
*2.56 (-5) 


1.99 (0) 
1.96 (0) 
1.94 (0) 
1.87 (0) 

1.76 (0) 
1.55 (0) 

1.07 (0) 
*1.07 (0) 

5.85 (-1) 
*5.84 (-1) 

1.77 (-1) 
*1.77 (-1) 

*5.70 (-3) 
*2.41 (-5) 


1.72 (0) 
1.70 (0) 
1.68 (0) 
1.61 (0) 
1.51 (0) 
1.32 (0) 

8.96 (-1) 
*8.94 (-1) 

4.73 (-1) 
*4.73 (-1) 

1.37 (-1) 
*1.38 (-1) 

*4.21 (-3) 
*1.78 (-5) 


1.51 (0) 
1.49 (0) 
1.47 (0) 
1.41 (0) 
1.31 (0) 
1.14 (0) 

7.56 (-1) 
*7.55 (-1) 

3.87 (-1) 
*3.87 (-1) 


1.30 (0) 
1.26 (0) 
1.20 (0) 
1.20 (0) 
1.12 (0) 
9.62 (-1) 

6.18 (-1) 
•6.17 (-1) 

3.02 (-1) 
*3.04 (-1) 


1.01 (0) 
9.95 (-1) 
9.78 (-1) 
9.29 (-1) 
8.53 (-1) 
7.20 (-1) 

4.37 (-1) 
*4.37 (-1) 

1.95 (-1) 
*1.95 (-1) 


6.66 (-1) 
6.52 (-1) 

6.38 (-1) 
5.99 (-1) 

5.39 (-1) 
4.37 (-1) 

2.35 (-1) 
*2.42 (-1) 

8.65 (-2) 
*8.65 (-2) 


1.23 (-1) 
1.17 (-1) 
1.11 (-D 
9.57 (-2) 
7.45 (-2) 
4.52 (-2) 

1.02 (-2) 
*9.31 (-3) 


3.13 (-2) 
2.83 (-2) 
2.56 (-2) 
1.90 (-2) 
1.15 (-2) 
4.24 (-3) 




*3.96 (-4) 




*1.06 (-3) 


*2.27 (-6) 




*1.06 (-1) 

*2.95 (-3) 
*1.14 (-5) 


*7.67 (-2) 

*1.79 (-3) 
*5.57 (-6) 


*4.22 (-2) 

*6.51 (-4) 
•1.13 (-6) 


*1.29 (-2) 

*7.40 r-5) 
*2.66 (-8) 


*6.65 (-6). 




















ft=2.0 


\ 
\ P 

\ 

r \ 

\ 





0.2 


0.5 


0.8 


1.0 


1.2 


1.5 


2.0 


5.0 


10.0 






0.01 

0.02 

0.05 

0.1 

0.2 

0.5 

1.0 

2.0 

5.0 
10.0 


7.01 (-1) 
6.94 (-1) 
6.86 (-1) 
6.65 (-1) 
6.31 (-1) 
5.67 (-1) 

4.13 (-1) 
*4. 14 (-1) 

2.44 (-1) 
*2.43 (-1) 

8.47 (-2) 
*8.44 (-2) 

*3.59 (-3) 
*1.89 (-5) 


6.96 (-1) 
6.89 (-1) 
6.81 (-1) 
6.60 (-1) 
6.26 (-1) 
5.63 (-1) 

4.09 (-1) 
*4. 10 (-1) 

2.41 (-1) 
*2.41 (-1) 

8.36 (-2) 
*8.33 (-2) 

*3.52 (-3) 
*1. 84 (-5) 


6.71 (-1) 
6.63 (-1) 
6.56 (-1) 
6.35 (-1) 
6.02 (-1) 
5.40 (-1) 

3.91 (-1) 
*3.90 (-1) 

2.28 (-1) 
*2.28 (-1) 

7.79 (-2) 
*7.79 (-2) 

*3. 17 (-3) 
♦1.61 (-5) 


6.27 (-1) 
6.20 (-1) 
6.14 (-1) 
5.93 (-1) 
5.61 (-1) 
5.02 (-1) 

3.60 (-1) 
*3.59 (-1) 

2.06 (-1) 
*2.07 (-1) 

6.85 (-2) 
*6.85 (-2) 

*2.60 (-3) 
*1.24 (-5) 


5.91 (-1) 
5.85 (-1) 
5.78 (-1) 
5.58 (-1) 
5.27 (-1) 
4.70 (-1) 

3.34 (-1) 
*3.34 (-1) 

1.89 (-1) 
*1.89 (-1) 

6.09 (-2) 
*6.09 (-2) 

•2. 16 (-3) 
*9. 55 (-6) 


5.52 (-1) 
5.45 (-1) 
5.39 (-1) 
5.20 (-1) 
4.90 (-1) 
4.35 (-1) 

3.05 (-1) 
*3.07 (-1) 

1.70 (-1) 
*1.70 (-1) 

5.29 (-2) 
*5.29 (-2) 

*1.72 (-3) 
*6.86 (-6) 


4.90 (-1) 
4.84 (-1) 
4.78 (-1) 
4.60 (-1) 
4.32 (-1) 
3.81 (-1) 

2.62 (-1) 
•2. 62 (-1) 

1.41 (-1) 
*1.40 (-1) 

4. 12 (-2) 
*4. 12 (-2) 

•1. 14 (-3) 
*3. 45 (-6) 


3.92 (-1) 
3.86 (-1) 
3.81 (-1) 
3.65 (-1) 
3.40 (-1) 
2.96 (-1) 

1.95 (-1) 
*1.95 (-1) 

9.74 (-2) 
*9.74 (-2) 


1.10 (-1) 
1.07 (-1) 
1.04 (-1) 
9.60 (-2) 
8.40 (-2) 
6.44 (-2) 

2.91 (-2) 
*3.00 (-2) 

7.83 (-3) 
*7.69 (-3) 


3.03 (-2) 
2.88 (-2) 
2.74 (-2) 
2.35 (-2) 
1.83 (-2) 
1.10 (-2) 

2.41 (-3) 
*2. 18 (-3) 




*2.35 (-4) 




*2.50 (-2) 

*4.77 (-4) 
*9.02 (-7) 


*6.16 (-4) 
*2.28 (-7) 
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Table 3. Unscattered component (lTr!a)I° (r,p,h) of the radiation flux calculated (a) using table 1 of this work and table 2 of 
reference [3] (numbers preceded by asterisks) and (b) using table 2 of this work (no asterisks) — Continued 













(iirla) 1° 


(p, K r) 












h=5.0 


\ 
\ P 
\ 

r \ 





0.2 


0.5 


0.8 


1.0 


1.2 


1.5 


2.0 


5.0 


10.0 






0.01 

0.02 

0.05 

0.1 

0.2 

0.5 

1.0 

2.0 

5.0 
10.0 


1.23 (-1) 
1.22 (-1) 
1.21 (-1) 
1.17 (-1) 
1.11 (-D 
1.01 (-1) 

7.44 (-2) 
*7.24 (-2) 

4.49 (-2) 
*4. 50 (-2) 

1.64 (-2) 
*1.66 (-2) 

*7.91 (-4) 
*5.09 (-6) 


1.23 (-1) 
1.22 (-1) 
1.21 (-1) 
1.17 (-1) 
1.11 (-D 
1.01 (-1) 

7.43 (-2) 
*7.24 (-2) 

4.48 (-2) 
*4.49 (-2) 

1.63 (-2) 
*1.66 (-2) 

*7.87 (-4) 
*5.05 (-6) 


1.22 (-1) 
1.21 (-1) 
1.20 (-1) 
1.16 (-1) 
1.10 (-1) 
9.96 (-2) 

7.35 (-2) 
*7.24 (-2) 

4.43 (-2) 
*4.44 (-2) 

1.60 (-2) 
*1.62 (-2) 

*7.66 (-4) 
*4.83 (-6) 


1.20 (-1) 
1.19 (-1) 
1.18 (-1) 
1.14 (-1) 
1.09 (-1) 
9.81 (-2) 

7.22 (-2) 
*7.20 (-2) 

4.33 (-2) 
*4.33 (-2) 

1.56 (-2) 
*1.57 (-2) 

*7.29 (-4) 
*4.45 (-6) 


1.19 (-1) 
1.18 (-1) 
1.16 (-1) 
1.13 (-1) 
1.07 (-1) 
9.66 (-2) 

7.10 (-2) 
*7. 14 (-2) 

4.25 (-2) 
*4.25 (—2) 

1.52 (-2) 
*1.52 (-2) 

*6.97 (-4) 
*4. 13 (-6) 


1.17 (-1) 
1.16 (-1) 
1.14 (-1) 
1.11 (-1) 
1.05 (-1) 
9.50 (-2) 

6.96 (-2) 
*7.03 (-2) 

4.14 (-2) 
*4. 14 (-2) 

1.47 (-2) 
*1.47 (-2) 

*6.60 (-4) 
*3. 78 (-6) 


1.14 (-1) 
1.12 (-1) 
1.11 (-D 
1.08 (-1) 
1.02 (-1) 
9.20 (-2) 

6.71 (-2) 
*6.80 (-2) 

3.97 (-2) 
*3.96 (-2) 

1.39 (-2) 
*1.38 (-2) 

*5.97 (-4) 
*3.20 (-6) 


1.07 (-1) 
1.06 (-1) 
1.05 (-1) 
1.01 (-1) 
9.60 (-2) 
8.62 (-2) 

6.23 (-2) 
*6.24 (-2) 

3.62 (-2) 
*3.62 (-2) 

1.23 (-2) 
*1.22 (-2) 

*4.83 (-4) 
*2.27 (-6) 


6.28 (-2) 
6.27 (-2) 
6.11 (-2) 
5.85 (-2) 
5.46 (-2) 
4.74 (-2) 

3.10 (-2) 
*3.14 (-2) 

1.53 (-2) 
*1.53 (-2) 

3.76 (-3) 
*3.77 (-3) 

*5. 71 (-5) 
*5.99 (-8) 


2.52 (-2) 
2.51 (-2) 
2.41 (-2) 

2.25 (-2) 
2.02 (-2) 
1.61 (-2) 

8.26 (-3) 
*8.50 (-3) 

2.71 (-3) 
*2.72 (-3) 




*2.94 (-4) 
*3.77 (-7) 







h = V..O 
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\ 
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1.5 
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10.0 






0.01 

0.02 

0.05 

0.1 

0.2 


3.13 (-2) 
3.09 (-2) 
3.06 (-2) 
2.97 (-2) 
2.83 (-2) 
2.56 (-2) 


3.12 (-2) 
3.09 (-2) 
3.06 (-2) 
2.97 (-2) 
2.83 (-2) 
2.56 (-2) 


3.12 (-2) 
3.09 (-2) 
3.06 (-2) 
2.97 (-2) 
2.82 (-2) 
2.55 (-2) 


3.11 (-2) 
3.08 (-2) 
3.04 (-2) 
2.95 (-2) 
2.81 (-2) 
2.54 (-2) 


3.10 (-2) 
3.06 (-2) 
3.03 (-2) 
2.94 (-2) 
2.80 (-2) 
2.53 (-2) 


3.08 (-2) 
3.05 (-2) 
3.02 (-2) 
2.93 (-2) 
2.79 (-2) 
2.52 (-2) 


3.06 (-2) 
3.03 (-2) 
3.00 (-2) 
2.91 (-2) 
2.76 (-2) 
2.50 (-2) 


3.01 (-2) 
2.98 (-2) 
2.95 (-2) 
2.86 (-2) 
2.72 (-2) 
2.45 (-2) 


2.51 (-2) 
2.48 (-2) 
2.45 (-2) 
2.37 (-2) 
2.24 (-2) 
2.00 (-2) 


1.57 (-2) 
1.55 (-2) 
1.53 (-2) 
1.46 (-2) 
1.36 (-2) 
1.18 (-2) 




0.5 


1.89 (-2) 
*1. 76 (-2) 


1.89 (-2) 
•1.76 (-2) 


1.89 (-2) 
•1.77 (-2) 


1.88 (-2) 
•1.79 (-2) 


1.87 (-2) 
*1.80 (-2) 


1.86 (-2) 
*1.81 (-2) 


1.84 (-2) 
*1.82 (-2) 


1.80 (-2) 
•1.83 (-2) 


1.43 (-2) 
•1.40 (-2) 


7.75 (-3) 
*7.86 (-3) 




1.0 


1. 15 (-2) 
*1. 15 (-2) 


1.15 (-2) 
•1.15 (-2) 


1.14 (-2) 

n. is (-2) 


1.14 (-2) 
•1. 14 (-2) 


1.13 (-2) 
*1.13 (-2) 


1.12 (-2) 
*1. 13 (-2) 


1.11 (-2) 
•1.11 (-2) 


1.08 (-2) 
•1.08 (-2) 


8.20 (-3) 
*8.20 (-3) 


3.82 (-3) 
*3.80 (-3) 




2.0 


4.21 (-3) 
*4. 33 (-3) 


4.21 (-3) 
*4.32 (-3) 


4.19 (-3) 
*4.30 (-3) 


4.16 (-3) 
*4.25 (-3) 


4.13 (-3) 
•4.21 (-3) 


4.09 (-3) 
•4.16 (-3) 


4.03 (-3) 
*4.08 (-3) 


3.90 (-3) 
•3.91 (-3) 


2.68 (-3) 
*2.67 (-3) 


9.31 (-4) 
*9.33 (-4) 




5.0 
10.0 


*2.08 (-4) 
*1.39 (-6) 


*2.08 (-4) 
•1.38 (-6) 


*2.06 (-4) 
•1.37 (-6) 


*2.04 (-4) 
•1.34 (-6) 


*2.04 (-4) 
•1.31 (-6) 


•2.01 (-4) 
*1.28 (-6) 


•1.93 (-4) 
•1.22 ( -6) 


*1.82 (-4) 
*1.10 (-6) 


*9.37 (-5) 
*3.53 (-7) 


*1.36 (-5) 
•1.26 (-8) 



Table 4. Comparison of convergence in calculating the un- 
scattered component (4ir/<r)I (T.p t h) in the case T=nt = l, for 
the three situations (A), (£), and (C) indicated in figure 8. 
For each situation the first column contains partial su?ns 



X! (l + i)gi(r)pi( Pl h) 

1 = 



and the second column contains sums 
S 9n(p,h)r\ 

71 = 





(A): P = 


0.5, fc=0.5 


(B): p- 


= l,h = l 


(C): P = 


= 2, h = l 


L,N 


Legendre 


r n series 


Legendre 


r n series 


Legendre 


r n -series 




sum 


sum 


sum 


sum 


sum 


sum 





0. 34031 


4. 5353 


0. 12316 


1.5118 


0.035740 


0. 66582 


1 


. 84669 


-2. 6711 


. 31658 


-0. 6170 


. 072679 


-0. 75149 


2 


. 99567 


3. 6122 


. 38273 


0. 9538 


. 068353 


0. 81932 


3 


. 99211 


-0. 4003 


. 37916 


.1441 


. 071785 


-0. 38701 


4 


1. 01738 


1. 6941 


. 38475 


.4713 


. 080105 


0. 33294 


5 


1. 02484 


0. 7540 


. 38774 


.3612 


. 079609 


-0. 02187 


6 


1. 02365 


1.1264 


. 38665 


.3932 


. 082707 


0. 12793 


7 


1. 02628 


0. 9943 


. 38793 


.3850 


. 084872 


. 07242 


8 


1. 02590 


1. 0367 


. 38753 


.3869 


. 084935 


. 09079 


9 


1.02615 


1. 0243 


. 38718 


.3865 


. 086201 


. 08529 


10 


1.02619 


1. 0276 


. 38727 


. 3865 


. 086348 


. 08679 


11 


1. 02655 


1. 0268 


. 38688 


.3865 


. 086408 


. 08641 


12 


1. 02677 


1. 0270 


. 38685 


.3865 


. 086602 


. 08650 


13 


1. 02685 


1. 0269 


. 38672 


.3865 


. 086496 


. 08648 



A physical interpretation of the rather wild 
oscillations in the first few terms of the r n series is 
that the zeroth term would be the result for a bare 
disk source (r = 0), from which the radiation ab- 
sorbed in the barrier of thickness r must be sub- 
tracted. This subtraction is accomplished by a 
succession of undershoots and overshoots which are 
quickly damped out if r is small enough. 



3.4. A«0 # P >1: 



Power Series Solution In or Near 
Plane of Disk 



If both source and detector are embedded in a 
homogeneous attenuating medium, the series (34) is 
valid except for h=0, where the concept of "barrier 
thickness" t becomes meaningless. As h->0, q n ->°° 
as l/h n , but at the same time r n ^0 as h n , while 
qo(pfi) is a finite but indeterminate quantity of the 
form { oo — oo } . 

The series (34) can still be used if the indeter- 
minacy is removed from the product q n r n by multiply- 
ing and dividing q n and r n , respectively, by h n y 
resulting in a series of the form 



/'(T,p,A) = (<r/4ir)2 2'»(p,A)(/i-l)», 

n=Q 



(47) 
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in which q' n (p,h) = h n q n (p,h) and p-l = T/h = pt/h= 
nh/h= the radius of the disk source measured in 
mean-free-paths of the attenuating medium, and 

For h=0, p^>l the expression in (37) reduces to 



go(p,0)=ln 



*> 



I Is 



4. Buildup Factor 



When both source and detector are embedded in 
a homogeneous medium, the secondary, or scattered 
radiation contribution I s is often taken into account 
by means of a buildup factor [6] 

B(»r) = {I +I s )II\ 

Such a buildup factor may then be used to evaluate 
the total flux / from a point isotropic source ac- 
cording to 

I=P+P=PB(jtr) = (c/4ar)B{ia') exp(-jur)/r 2 . (49) 
In references [8] and [9] B (jur) is expressed as a set 

• SEMI-INFINITE MEDIUM 

ISOTROPIC EXTENDED SOURCE 



of polynomial coefficients fi i} in the notation of 
reference [9], as defined in equation (2), with N=3. 
Now, since the integration of (33) was performed 
term by term over the expansion 

exp (-t/cos 0)=exp (-fir)=^(-l) n ( f xr) n /n\ ) (50) 

no new integrands appear if these series terms are 
multiplied by terms fi i (fir) i of the polynomial repre- 
sentation of B(fxr). if this multiplication is per- 
formed, the total radiation flux as defined in (49) 
may be written as the series 



71 = 



(51) 



in which r=nh. The modified buildup coefficients 
fi' n are generated from /3/s, such as are given in 
[8] or [9], as 



^=^(-iy{n\/(n-i)l}^ 



(52) 



in which n\/(n— i)!=0 for i^>n. For N=S the 
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Figure 5. Some circular aperture or disk source situations to which results are applicable. 
(a): I 0j rl a calculable using pi(p,h), (Ji(t) coefficients, (b) : I°-{-I a calculable using Pj,qj(p,h) 
coefficients, (a-f): 1° calculable using either pi(p,h) and gi°, or qj(p,h) coefficients. 
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first five modified buildup coefficients are 
flS=i8b— 2A+2 - 1A 



and 



/3£=A>-3ft+3- 2&-3.2 -103 



0;=/3 o -4ft+4 • 30 2 -4 • 3 • 20 3 . 



(53) 



Such a set of 0^ coefficients, once generated, 
may be applied as in (51) to either the disk-source 
q n (P)hy& in table 3 of this work or the rectangular 
source g w (a,&)'s in table 3 of reference [3]. 

5. Discussion 

Situations to which results in this paper may be 
applied are schematized in figure 5. 

The Legendre method of section 2 can be used 
to calculate I=I°-\-I s in such cases as figure 5a 
in which the source is distributed in a uniform layer 
(or layers) parallel to the exit surface of the medium. 
If gi angular harmonic coefficients are used which 
have been calculated assuming an infinite medium, 
the result will be an overestimate because of the 
depression, in the realistic case, of flux near a 
boundary. This flux depression results from the 
one-way nature of the radiation traversal across 
the boundary. For gamma rays the infinite medium 
g/s may be a rather good approximation, whereas 
for neutrons more realistic gjs including the boundary 
effect may have to be used. 

The unscattered component 1° may in principle 
be calculated for all situations in figure 5a-f using 
either the Legendre pi(p,h)'s of section 2 or the 
q n (p,hys of section 3, but in practice the Legendre 
series (6) converges better for barriers of the order 
of one mean free path or thicker while the q n (p,h) 
series (34) converges better for barriers of less than 
a mean free path. 

Adaptation of polynomial representations of point 
isotropic source buildup factors, as suggested in 
section 4 for use with the q n (p,hYs of section 3 for 
computation of I=I 0J rP, is strictly valid only 
for the situation in figure 5b, although some qualita- 
tive results may be obtained in the other situations 
depending on how little the physical situation 
departs from figure 5b. 
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